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ABSTRACT 

We  describe  a  standard  star  catalog  constructed  using  multiple  SDSS  photometric  observations  (at  least  four  per 
band,  with  a  median  of  10)  in  the  ugriz  system.  The  catalog  includes  1 .01  million  nonvariable  unresolved  objects  from 
the  equatorial  stripe  82  (|6j2ooo.o|  <  1.266°)  in  the  right  ascension  range  20h34ra-4h00m  and  with  the  corresponding 
r-band  (approximately  Johnson  K-band)  magnitudes  in  the  range  14-22.  The  distributions  of  measurements  for  in¬ 
dividual  sources  demonstrate  that  the  photometric  pipeline  correctly  estimates  random  photometric  errors,  which  are 
below  0.01  mag  for  stars  brighter  than  19.5, 20.5, 20.5,  20,  and  18.5  in  ugriz,  respectively  (about  twice  as  good  as  for 
individual  SDSS  runs).  Several  independent  tests  of  the  internal  consistency  suggest  that  the  spatial  variation  of  pho¬ 
tometric  zero  points  is  not  larger  than  ^0.01  mag  (rms).  In  addition  to  being  the  largest  available  data  set  with  optical 
photometry  internally  consistent  at  the  ~1%  level,  this  catalog  provides  a  practical  definition  of  the  SDSS  photomet¬ 
ric  system.  Using  this  catalog,  we  show  that  photometric  zero  points  for  SDSS  observing  runs  can  be  calibrated  within 
a  nominal  uncertainty  of  2%  even  for  data  obtained  through  1  mag  thick  clouds,  and  we  demonstrate  the  existence  of 
He  and  H  white  dwarf  sequences  using  photometric  data  alone.  Based  on  the  properties  of  this  catalog,  we  conclude 
that  upcoming  large-scale  optical  surveys  such  as  the  Large  Synoptic  Survey  Telescope  will  be  capable  of  delivering 
robust  1%  photometry  for  billions  of  sources. 

Key  words:  catalogs  —  instrumentation:  photometers  —  methods:  data  analysis  —  standards  —  surveys  — 
techniques:  photometric 
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1.  INTRODUCTION 

Astronomical  optical  photometric  data  are  usually  calibrated 
using  sets  of  standard  stars  whose  brightness  is  known  from 
previous  work.  The  most  notable  modem  optical  standard  star 
catalogs  are  the  Landolt  standards  (Landolt  1992)  and  Stetson 
standards  (  Stetson  2000, 2005 ).  Both  are  reported  on  the  Johnson- 
Kron-Cousins  system  (Landolt  1983  and  references  therein). 
The  Landolt  catalog  provides  magnitudes  accurate  to  L  I  %  in 
the  UB  VR1  bands  for  ~500  stars  in  the  V magnitude  range  11.5— 
1 6.  Stetson  has  extended  Landolt’s  work  to  fainter  magnitudes  and 
provided  the  community  with  ~  1  %— 2%  accurate  magnitudes  in 
the  RFR/ bands  for  ~  15, 000  starsinthe  magnitude  range  V  51  20. 
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Most  stars  from  both  sets  are  distributed  along  the  celestial  equa¬ 
tor,  which  facilitates  their  use  from  both  hemispheres. 

The  data  obtained  by  the  Sloan  Digital  Sky  Survey  (SDSS; 
York  et  al.  2000)  can  be  used  to  extend  the  work  by  Landolt  and 
Stetson  to  even  fainter  levels  and  to  increase  the  number  of  stan¬ 
dard  stars  to  over  1  million.  In  addition,  SDSS  has  designed  its 
own  photometric  system  (ugriz;  Fukugita  et  al.  1996),  which  is 
now  in  use  at  many  observatories  worldwide.  This  widespread  use 
of  the  ugriz  photometric  system  motivates  the  construction  of  a 
large  standard  star  catalog  with  ~  1%  accuracy.  As  part  of  its  im¬ 
aging  survey,  SDSS  has  obtained  many  scans  in  the  so-called 
stripe  82  region,  which  is  defined  by  |<5j2ooo.o|  <  1-266°  and  right 
ascension  approximately  in  the  range  20h-4h.  These  repeated 
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observations  can  be  averaged  to  produce  more  accurate  pho¬ 
tometry  than  the  nominal  2%  single-scan  accuracy  (Ivezic  et  al. 
2004a). 

The  catalog  and  methods  presented  here  have  some  similarity 
to  an  effort  by  Padmanabhan  et  al.  (2007),  who  developed  a  new 
calibration  algorithm  that  simultaneously  solves  for  the  calibra¬ 
tion  parameters  and  relative  stellar  fluxes  using  overlapping  SDSS 
observations  (the  so-called  tibercalibration  method).  The  algo¬ 
rithm  decouples  the  problem  of  “relative”  calibrations  (i.e.,  pro¬ 
ducing  an  internally  consistent  system)  from  that  of  “absolute” 
calibrations  (i.e.,  typing  the  internal  system  to  a  physical  flux 
scale):  the  absolute  calibration  is  reduced  to  determining  a  few 
numbers  for  the  entire  survey.  Here  we  also  decouple  relative  and 
absolute  calibrations  and  use  overlapping  observations.  The  main 
difference  between  their  work  and  this  paper  is  that  they  are  con¬ 
cerned  about  calibrating  the  entire  SDSS  survey  (^8500  deg2;  most 
of  the  surveyed  area  has  at  most  two  overlapping  observations), 
while  we  concentrate  here  on  a  much  smaller  area  (~300  deg2) 
with  an  average  of  10  overlapping  observations.  We  also  deter¬ 
mine  flat-field  corrections  (relative  to  the  “standard”  survey  re¬ 
ductions)  using  different  methods:  Padmanabhan  et  al.  minimize 
errors  in  relative  photometry  of  multiply  observed  stars,  while 
we  require  that  the  stellar  locus  remains  fixed  in  multidimensional 
color  space.  An  advantage  of  the  catalog  presented  here  is  better 
averaging  of  various  photometric  errors  thanks  to  a  larger  num¬ 
ber  of  observations,  which  comes  at  the  expense  of  a  much  smaller 
cataloged  area.  It  is  encouraging  that  the  results  of  these  two  com¬ 
plementary  approaches  agree  in  the  regions  of  sky  common  to  both 
catalogs  at  the  claimed  level  of  accuracy  (~1%). 

Additional  motivation  for  the  analysis  of  repeated  scans  and 
their  impact  on  photometric  accuracy  comes  from  upcoming  large- 
scale  optical  surveys  such  as  the  Dark  Energy  Survey  ( Flaugher 
et  al.  2006),  Pan-STARRS  (Kaiser  et  al.  2002),  and  the  Large 
Synoptic  Survey  Telescope  (LSST;  Tyson  2002).  For  example, 
the  LSST  science  requirements  document19  calls  for  a  photometric 
system  that  is  internally  consistent  across  the  sky  at  the  1%  level. 
The  SDSS  stripe  82  repeated  scans  can  be  used  to  gauge  the  plaus¬ 
ibility  of  delivering  such  a  system. 

We  describe  the  construction  and  testing  of  a  standard  star  cat¬ 
alog  in  §  2  and  illustrate  its  many  uses  in  §  3.  We  discuss  our  re¬ 
sults  in  §  4. 

2.  THE  CONSTRUCTION  OF  THE  SDSS  STRIPE  82 
STANDARD  STAR  CATALOG 

2.1.  Overview  of  SDSS  Imaging  Data 

SDSS  uses  a  dedicated  2.5  m  telescope  (Gunn  et  al.  2006)  to 
provide  homogeneous  and  deep  (r  <  22.5)  photometry  in  five 
bandpasses  (Fukugita  et  al.  1996;  Gunn  et  al.  1998, 2006;  Smith 
et  al.  2002;  Hogg  et  al.  2001)  repeatable  to  0.02  mag  (rms  for 
sources  not  limited  by  photon  statistics;  Ivezic  et  al.  2003)  and 
with  a  zero-point  uncertainty  of  ~0.02-0.03  (Ivezic  etal.  2004a). 
The  survey  sky  coverage  of  close  to  10,000  deg2  in  the  northern 
Galactic  cap  and  ~300  deg2  in  the  southern  Galactic  hemisphere 
will  result  in  photometric  measurements  for  well  over  1 00  million 
stars  and  a  similar  number  of  galaxies.20  Astrometric  positions  are 
accurate  to  better  than  0.1"  per  coordinate  (rms)  for  sources  with 
r  <  20.5  mag  (Pier  et  al.  2003),  and  the  morphological  infonna- 


19  Available  at  http://www.lsst.org/Science/lsst_baseline.shtml. 

20  The  recent  Data  Release  5  (J.  K.  Adelman-McCarthy  et  al.  2007,  in  prepa¬ 
ration)  lists  photometric  data  for  2 15  million  unique  objects  observed  in  8000  deg2 
of  sky;  see  http://www.sdss.org/dr5/. 


tion  from  the  images  allows  reliable  star-galaxy  separation  to  r  ~ 
21.5  mag  (Lupton  et  al.  2002;  Scranton  et  al.  2002). 

Data  from  the  imaging  camera  (30  photometric,  12  astromet¬ 
ric,  and  2  focus  CCDs;  Gunn  et  al.  1998)  are  collected  in  drift- 
scan  mode.  The  images  that  correspond  to  the  same  sky  location 
in  each  of  the  five  photometric  bandpasses  (these  five  images  are 
collected  over  ~5  minutes,  with  54  s  for  each  exposure)  are 
grouped  together  for  simultaneous  processing  as  a  field.  A  field 
is  defined  as  a  36  s  (1361  pixels,  or  9';  see  Stoughton  et  al.  2002) 
long  and  2048  pixels  wide  (13')  stretch  of  drift-scanning  data  from 
a  single  column  of  CCDs  (sometimes  called  a  “scan  line”;  for 
more  details,  see  Stoughton  et  al.  2002;  Abazajian  et  al.  2003, 
2004, 2005;  Adelman-McCarthy  et  al.  2006).  Each  of  the  six  scan 
lines  (together  called  a  “strip”)  is  13'  wide.  The  12  interleaved 
scan  lines  (or  two  strips)  are  called  a  “stripe”  (~2.5°  wide). 

2.2.  The  Photometric  Calibration  of  SDSS  Imaging  Data 

SDSS  2.5  m  imaging  data  are  photometrically  calibrated  using  a 
network  of  calibration  stars  obtained  in  1520  41 .5  x  41 .5  arcmin2 
transfer  fields,  called  “secondary  patches.”  These  patches  are  po¬ 
sitioned  throughout  the  survey  area  and  are  calibrated  using  a  pri¬ 
mary  standard  star  network  of  158  stars  distributed  around  the 
northern  sky  (so-called  USNO  standards;  Smith  et  al.  2002).  The 
primary  standard  star  network  is  tied  to  an  absolute  flux  system  by 
the  single  F0  subdwarf  starBD  +17  4708,  whose  absolute  fluxes 
in  the  SDSS  filters  are  taken  from  Fukugita  et  al.  (1996).  The  sec¬ 
ondary  patches  are  grouped  into  sets  of  four  and  are  observed  by 
the  Photometric  Telescope  (PT;  Tucker  et  al.  2006)  in  parallel 
with  observations  of  the  primary  standards.  A  set  of  four  patches 
spans  all  12  scan  lines  of  a  survey  stripe  along  the  width  of  the 
stripe,  and  the  sets  are  spaced  along  the  length  of  a  stripe  at 
roughly  15°  intervals,  which  corresponds  to  1  hr  of  scanning  at 
the  sidereal  rate. 

SDSS  2.5  m  magnitudes  are  reported  on  the  “natural  system” 
of  the  2.5  m  telescope  defined  by  the  photon-weighted  effective 
wavelengths  of  each  combination  of  SDSS  filter,  CCD  response, 
telescope  transmission,  and  atmospheric  transmission  at  a  refer¬ 
ence  air  mass  of  1.3  as  measured  at  APO.21  The  magnitudes  are 
referred  to  as  the  ugriz  system  (which  differs  from  the  “primed” 
system,  u'g'r'i'z' ,  that  is  defined  by  the  USNO  standards).22  The 
reported  magnitudes23  are  corrected  for  the  atmospheric  extinction 
(using  simultaneous  observations  of  standard  stars  by  the  PT )  and 
thus  correspond  to  measurements  at  the  top  of  the  atmosphere24 
(except  for  the  fact  that  the  atmosphere  has  an  impact  on  the  wave¬ 
length  dependence  of  the  photometric  system  response).  The  mag¬ 
nitudes  are  reported  on  the  AB  system  (Oke  &  Gunn  1983), 
defined  such  that  an  object  with  a  specific  flux  of  Fv  =  3631  Jy 
has  m  =  0  (i.e.,  an  object  with  F„=  const,  has  an  AB  mag¬ 
nitude  equal  to  the  Johnson  V magnitude  at  all  wavelengths).  In 
summary,  given  a  specific  flux  of  an  object  at  the  top  of  the 
atmosphere,  Fj,(/i),  the  reported  SDSS  2.5  m  magnitude  in  a  given 


21  Transmission  curves  for  the  SDSS  2.5  m  photometric  system  are  available 
at  http://www.sdss.org  /dr5/instmments/imager. 

22  For  the  subtle  effects  that  led  to  this  distinction,  see  Stoughton  et  al.  (2002), 
Smith  et  al.  (2002),  and  http://www.sdss.org/dr5/algorithms/fluxcal.html. 

23  SDSS  uses  a  modified  magnitude  system  ( Lupton  et  al.  1 999),  which  is  vir¬ 
tually  identical  to  the  standard  astronomical  Pogson  magnitude  system  at  high 
signal-to-noise  ratios  relevant  here. 

24  The  same  atmospheric  extinction  correction  is  applied  irrespective  of  the 
source  color;  within  the  relevant  air-mass  range  (X  <  1 .5),  the  systematic  errors 
this  introduces  are  less  than  1%  for  stars  with  effective  temperatures  in  the  3000- 
20,000  K  range.  For  objects  with  significantly  different  spectral  energy  distribu¬ 
tions,  such  as  quasars  and  supemovae,  the  errors  may  be  several  times  larger.  See 
also  Maiz  Apellaniz  (2006). 
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band,  b  =  (u,  g ,  r,  i,z),  corresponds  to  (modulo  random  and  sys¬ 
tematic  errors,  which  are  discussed  later) 


m  =  —2.5  login  f - - — V  (1) 

SI°  V3631  Jy  /  W 

where 

Fb  =  J  Fv(X)(f>b(X)  dk.  (2) 


Here  <f>b(F)  is  the  normalized  system  response^5  for  the  given 
band, 


4>h(F) 


rlSbU) 

J  A~lSh(A)  dX  ’ 


(3) 


with  the  overall  atmosphere  plus  system  throughput,  5*  (A),  avail¬ 
able  from  the  Web  site  given  in  footnote  21  for  the  SDSS 

system  is  shown  in  Fig.  6,  discussed  below;  see  also  §  2.5.2].  We 
reiterate  that  the  normalization  of  reported  magnitudes  corresponds 
to  a  source  at  the  top  of  the  atmosphere,  while  the  throughput  <fib(F) 
includes  the  transmission  of  a  standard  atmosphere  at  a  fiducial 
air  mass  of  1.3.  Note  also  that  it  is  only  the  shape  of  5*(A),  and 
not  its  overall  normalization,  that  needs  to  be  known  to  compute 
expected  SDSS  magnitudes  of  a  source  with  a  given  Fv(a).  That 
is,  the  SDSS  photometric  system  is  fully  defined  by  the  five  scale¬ 
less  functions  </>a(A)  (by  definition,  J  (j>bdX  =  l;  see  eq.  [3]).  In 
reality,  for  each  ugriz  band  there  are  six  devices  in  the  SDSS 
camera  (Gunn  et  al.  1998)  whose  fb  values  are  slightly  different 
(see  §  2.5.2). 

The  quality  of  SDSS  photometry  stands  out  among  the  avail¬ 
able  large-area  optical  sky  surveys  (Ivezic  et  al.  2003,  2004a; 
Sesar  et  al.  2006).  Nevertheless,  the  achieved  accuracy  is  occa¬ 
sionally  worse  than  the  nominal  0.02-0.03  mag  (rms  scatter  for 
sources  not  limited  by  photon  statistics).  Typical  causes  of  sub¬ 
standard  photometry  include  an  incorrectly  modeled  point-spread 
function  (PSF;  usually  due  to  fast  variations  of  atmospheric  see¬ 
ing  or  the  lack  of  a  sufficient  number  of  the  isolated  bright  stars 
needed  for  modeling  the  PSF),  unrecognized  changes  in  atmo¬ 
spheric  transparency,  errors  in  photometric  zero-point  calibration, 
effects  of  crowded  fields  at  low  Galactic  latitudes,  an  undersam¬ 
pled  PSF  in  excellent  seeing  conditions  (510.8";  the  pixel  size  is 
0.4"),  incorrect  flat-field  or  bias  vectors,  and  scattered  light  erro¬ 
neously  included  in  the  flat  field.  Such  effects  can  conspire  to  in¬ 
crease  the  photometric  errors  to  levels  as  high  as  0.05  mag  (with  a 
frequency,  at  that  error  level,  of  roughly  one  field  per  thousand). 
However,  when  multiple  scans  of  the  same  sky  region  are  avail¬ 
able,  many  of  these  errors  can  be  minimized  by  properly  averag¬ 
ing  the  photometric  measurements. 


2.3.  The  Choice  of  Cataloged  Magnitudes 

The  SDSS  photometric  pipeline  ( photo;  Lupton  et  al.  2002) 
measures  several  types  of  magnitudes,  including  aperture,  PSF, 
and  model  magnitudes.  Here  we  briefly  describe  each  type  of 
magnitude  (for  more  details,  see  Stoughton  et  al.  2002  and  the 
SDSS  Web  site)26  and  justify  the  choice  of  PSF  magnitudes  for 
catalog  construction. 


25  The  term  A-1  in  eq.  (3)  reflects  the  fact  that  CCDs  are  photon-counting  de¬ 
vices;  for  more  details,  see  Appendix  A  in  Maiz  Apellaniz  (2006 ). 

26  Available  at  http://www.sdss.org. 


2.3.1.  Aperture  Magnitudes 

Aperture  magnitudes  computed  by  photo  are  based  on  the  flux 
contained  within  the  aperture  with  a  radius  of  7.43".  While  an 
aperture  magnitude  is  the  most  robust  flux  estimate  for  unsatu¬ 
rated  sources  at  the  bright  end  (because  it  is  essentially  seeing- 
independent),  these  magnitudes  do  not  have  good  noise  proper¬ 
ties  at  the  faint  end  where  sky  noise  dominates  (e.g.,  for  a  given 
maximum  photometric  error,  PSF  magnitudes  reach  1-1.5  mag 
fainter  than  aperture  magnitudes).  In  order  to  improve  the  depth 
of  the  standard  star  catalog,  we  opt  not  to  use  aperture  magnitudes, 
except  for  quality  tests  at  the  bright  end. 

2.3.2.  Point-Spread  Function  Magnitudes 

The  PSF  flux  is  computed  using  the  PSF  as  a  weighting  func¬ 
tion.  While  this  flux  is  optimal  for  faint  point  sources  (in  partic¬ 
ular,  it  is  vastly  superior  to  aperture  photometry  at  the  faint  end), 
it  is  also  sensitive  to  inaccurate  PSF  modeling  as  a  function  of  po¬ 
sition  and  time.  Even  in  the  absence  of  atmospheric  variations,  the 
SDSS  telescope  and  camera  optics  deliver  images  whose  F  WHMs 
vary  by  up  to  1 5%  from  one  side  of  a  CCD  to  the  other;  the  worst 
effects  are  seen  in  the  chips  farthest  from  the  optical  axis.  More¬ 
over,  since  the  atmospheric  seeing  varies  with  time,  the  delivered 
image  quality  is  a  complex  two-dimensional  function  even  on 
the  scale  of  a  single  frame.  Without  accounting  for  this  spatial 
variation,  the  PSF  photometry  would  have  errors  up  to  0.10- 
0.15  mag.  The  description  of  the  PSF  is  also  critical  for  star-galaxy 
separation  and  for  unbiased  measures  of  the  shapes  of  nonstellar 
objects. 

The  SDSS  imaging  PSF  is  modeled  heuristically  in  each  band 
and  each  camera  column  using  a  Karhunen-Loeve  ( KL)  transform 
(Lupton  et  al.  2002).  Using  stars  brighter  than  roughly  20  mag, 
the  PSF  from  a  series  of  five  frames  is  expanded  into  eigenimages, 
and  the  first  three  terms  are  retained.  The  variation  of  these  coef¬ 
ficients  is  then  fit  up  to  a  second-order  polynomial  in  each  chip 
coordinate.  The  failure  of  this  KL  expansion,  typically  due  to  an 
insufficient  number  of  stars  or  an  exceedingly  complex  PSF,  re¬ 
sults  in  occasional  problems  with  PSF  photometry.  The  main 
failure  mode  is  inaccurate  determination  of  aperture  corrections 
which  statistically  tie  PSF  magnitudes  to  aperture  magnitudes 
using  bright  stars. 


2.3.3.  Model  Magnitudes 

Just  as  the  PSF  magnitudes  are  optimal  measures  of  the  fluxes 
of  stars,  the  optimal  measure  of  the  flux  of  a  galaxy  uses  a  matched 
galaxy  model.  With  this  in  mind,  the  photometric  pipeline  fits  two 
models  to  the  two-dimensional  image  of  each  object  in  each  band: 
a  pure  de  Vaucouleurs  profile  and  a  pure  exponential  profile.27  Be¬ 
cause  the  models  are  convolved  with  a  double-Gaussian  fit  to  the 
PSF,  the  seeing  effects  are  accounted  for.  Aperture  corrections  are 
applied  to  make  these  model  magnitudes  equal  the  PSF  magni¬ 
tudes  in  the  case  of  an  unresolved  object. 

2.3.4.  The  Choice  of  Magnitudes  for  the  Standard  Star  Catalog 

A  comparison  between  aperture,  PSF,  and  model  magnitudes 
for  unresolved  sources  is  done  automatically  for  every  SDSS  ob¬ 
serving  run  ( runQA  pipeline;  Ivezic  et  al.  2004a).  An  analysis  of 
over  200  runs  indicates  that  model  magnitudes  are  more  robust 
than  PSF  magnitudes:  PSF  magnitudes  show  systematic  offsets 
of  0.05  mag  from  aperture  magnitudes  3  times  more  often  than 
model  magnitudes  (roughly  once  per  thousand  fields),  as  model 


27  For  more  details,  see  http://www.sdss.org/dr5/algorithms/photometry.html. 
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fits  have  more  degrees  offreedom.  On  the  other  hand,  an  analysis 
of  repeated  scans  indicates  that  estimates  of  photometric  errors 
by  the  photometric  pipeline  are  more  accurate  for  PSF  magnitudes 
(agreeing  at  the  10%  level  with  the  measured  values;  see  Ivezic 
et  al.  2003;  Scranton  et  al.  2005 )  than  for  model  magnitudes  (which 
are  smaller  than  the  measured  values  by,  typically,  30%-50%). 
Because  the  rejection  of  likely  variable  sources,  which  relies  on 
accurate  photometric  error  estimates,  is  an  important  step  in  the 
construction  of  the  standard  star  catalog  (see  below),  we  choose  to 
use  PSF  magnitudes  to  construct  the  catalog. 

2.4.  Catalog  Construction 

Using  58  SDSS-1  runs  from  stripe  82  (approximately  20h  < 
onooo.o  <  4h  and  |<5j2ooo.o|  <  1.266°,  but  not  all  runs  extend  over 
the  entire  right  ascension  range)  obtained  in  mostly  photometric 
conditions  (as  indicated  by  the  PT  calibration  residuals,  infrared 
cloud  camera,28  and  tests  performed  by  the  nmQA  pipeline),  can¬ 
didate  standard  stars  from  each  run  are  selected  by  requiring  the 
following: 

1 .  That  objects  are  classified  as  STAR  (based  on  the  difference 
between  model  and  PSF  magnitudes);  this  morphological  classifi¬ 
cation  really  means  unresolved  (point)  sources  (e.g.,  quasars  are 
also  included). 

2.  That  they  have  quoted  photometric  errors  in  the  PSF  mag¬ 
nitude  (as  computed  by  the  photometric  pipeline)  smaller  than 
0.05  mag  in  at  least  one  band. 

3.  That  the  processing  flags  BRIGHT,  SATUR,  BLENDED, 
or  EDGE  are  not  set  in  any  band.29 

These  criteria  select  unsaturated  point  sources  with  sufficiently 
high  signal-to-noise  ratio  per  single  observation  to  approach  final 
photometric  errors  of  0.02  mag  or  smaller. 

After  matching  all  detections  of  a  single  source  across  runs 
(using  a  1"  matching  radius),  various  photometric  statistics  such 
as  unweighted  mean,  median,  and  their  standard  errors;  rms  scat¬ 
ter;  number  of  observations;  and  y2  per  degree  of  freedom  are 
computed  for  magnitudes  in  each  band.  We  use  errors  reported  by 
the  photometric  pipeline  to  compute  \2  and  note  that  systematic 
errors  common  to  all  runs  do  not  contribute  to  its  value.  This  initial 
catalog  of  multiepoch  observations  includes  1 .4  million  point 
sources  with  at  least  four  observations  in  each  of  the  g,  r,  and  i 
bands.  The  median  number  of  observations  per  source  and  band 
is  10,  and  the  total  number  of  photometric  measurements  is 
~57  million. 

The  distributions  of  the  median  magnitudes,  their  standard  er¬ 
rors,  x2i  and  the  number  of  observations  for  a  subset  of  these 
sources  are  shown  in  Figure  1.  The  random  errors  in  the  median 
magnitude  [computed  as  0.928*1QR/(jV  —  1)1/2,  where  IQR  is 
the  25%-75%  interquartile  range  of  the  individual  measurement 
distribution  and  N  is  the  number  of  measurements;  note  that  the 
error  of  the  median  for  a  Gaussian  distribution  is  25%  larger  than 
the  error  of  the  mean;  Lupton  1993]  are  below  0.01  mag  at  the 
bright  end.  These  errors  are  reliably  computed  by  the  photomet¬ 
ric  pipeline,  as  indicated  by  the  y2  distributions.  The  distribu¬ 
tions  of  these  sources  in  color-magnitude  and  color-color  diagrams, 
constructed  using  median  magnitudes,  are  shown  in  Figure  2. 
For  a  detailed  interpretation  of  these  diagrams,  see  Lenz  et  al. 
(1998),  Fan  (1999),  Finlator  et  al.  (2000),  Helmi  et  al.  (2003), 


28  For  more  details  about  the  camera,  see  http://www.apo.nmsu.edu/Telescopes/ 
SDSS/eng.papers/19910801/19910801.html  and  Hogg  et  al.  (2001). 

29  For  more  details  about  processing  flags,  see  http://www.sdss.org/dr5/ 
products /catalogs /flags,  html  and  Stoughton  et  al.  (2002). 


and  Ivezic  et  al.  (2006 ).  It  is  evident  that  the  sample  is  dominated 
by  stars. 

2.4.1.  Selection  of  Candidate  Standard  Stars 

Adopted  candidate  standard  stars  must  have  at  least  four  obser¬ 
vations  and,  to  avoid  variable  sources,  %2  less  than  3  in  the  gri 
bands  (the  same  requirements  are  later  applied  in  the  u  and  z  bands 
when  using  the  catalog  for  calibration,  as  we  discuss  further  be¬ 
low).  We  also  limit  the  right  ascension  to  the  range  from  20h  34m 
to  4h  00m,  which  provides  a  simple  areal  definition  (together  with 
feooo.ol  <  1  -266°)  of  a  282  deg2  large  rectangular  region,  while 
excluding  only  a  negligible  fraction  of  stars.  With  the  final  con¬ 
dition  that  the  standard  error  for  the  mean  magnitude  in  the  r  band 
is  smaller  than  0.05  mag,  these  requirements  result  in  a  catalog 
with  slightly  over  1  million  sources  (1,006,849)  with  42  million 
photometric  measurements.  Of  those,  563,908  have  a  random  er¬ 
ror  for  the  median  magnitude  in  the  r  band  smaller  than  0.01  mag, 
and  for  405,954  stars  this  is  true  in  all  three  of  the  gri  bands.  Sub¬ 
sets  of  92,905  and  290,299  stars  satisfy  these  requirements  in  the 
ugri  and  griz  bands,  and  91,853  stars  satisfy  this  in  all  five  bands. 
The  distributions  of  candidate  standard  stars  that  satisfy  the  above 
selection  criteria  in  all  five  bands  in  color-magnitude  and  color- 
color  diagrams  are  shown  in  Figure  3. 

For  comparison,  the  distribution  of  sources  that  were  rejected 
as  variable  (\2  greater  than  3  in  at  least  one  of  the  gri  bands)  in 
color-magnitude  and  color-color  diagrams  is  shown  in  Figure  4. 
As  is  evident  from  a  comparison  with  Figure  3,  the  distribution 
of  variable  sources  in  the  color-color  diagrams  is  markedly  dif¬ 
ferent  from  that  of  nonvariable  sources.  It  is  especially  striking 
how  low-redshift  (z  <  2.2)  quasars  are  easily  detected  by  their 
variability  (for  more  details,  see  Ivezic  et  al.  2004b).  However,  it 
is  fairly  certain  that  not  all  variable  sources  are  recognized  as  such 
because  of  the  limited  number  of  repeated  observations  (~10). 
For  example,  an  eclipsing  binary  with  a  much  shorter  eclipse  du¬ 
ration  than  the  orbital  period  could  easily  escape  detection.  Anal¬ 
ysis  of  the  variable  subsample  is  presented  in  a  companion  paper 
(B.  Sesar  et  al.  2007,  in  preparation). 

The  sky  density  of  all  the  sources  and  those  selected  as  non¬ 
variable  are  shown  in  Figure  5.  At  high  Galactic  latitudes  (\b\  ~ 
60°)  the  fraction  of  candidate  standard  star  sources  is  ~80%. 

2.5.  Systematic  Photometric  Errors 

Photometric  errors  computed  by  the  photometric  pipeline  pro¬ 
vide  a  good  estimate  of  random  errors  in  SDSS  photometry,  as 
demonstrated  by  the  y2  distributions  shown  in  Figure  1.  How¬ 
ever,  the  measurements  are  also  subject  to  systematic  errors  such 
as  spatial  dependence  of  the  internal  zero  points  (calibration  er¬ 
rors)  and  the  overall  deviations  of  the  internal  SDSS  zero  points 
from  an  AB  magnitude  scale.  Formally,  the  true  AB  magnitude 
of  an  object  (defined  by  eq.  [1])  in  a  given  band,  wtrue,  can  be  ex¬ 
pressed  as 

^ true  ni cat-!"  ( h •  A . ,  decl . )  — |—  Am ,  (4) 

where  ?wcat  is  the  cataloged  magnitude,  fy,(R.A.,  decl.)  describes 
the  spatial  variation  of  the  internal  zero-point  error  around  A„, 
(thus,  the  average  of  bm  over  the  cataloged  area  is  zero  by  con¬ 
struction),  and  Am  is  the  overall  (spatially  independent)  devia¬ 
tion  of  the  internal  SDSS  system  from  a  perfect  AB  system  (the 
five  values  of  A„,  are  equal  for  all  the  cataloged  objects).  Here  we 
ignore  systematic  effects,  e.g.,  device  nonlinearity  and  bandpass 
variations  between  different  camera  columns,  which  depend  on 
individual  source  properties  such  as  brightness  and  color  (but  see 
§  2.5.2  below). 
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Fig.  1. — Median  magnitude  error  as  a  function  of  magnitude  (left  column,  symbols),  the  y2  per  degree  of  freedom  distribution  (middle  column,  solid  line),  and  the 
number  of  observations  in  each  band  (right  column)  for  candidate  standard  stars  from  SDSS  stripe  82  (selected  by  y2  <  3  in  the  gri  bands).  The  solid  lines  in  the  left 
column  show  cumulative  magnitude  distribution  normalized  to  1 0  at  the  faint  end.  The  dashed  lines  in  the  middle  column  show  the  X  2  Per  degree  of  freedom  distribution 
for  a  Gaussian  error  distribution  and  9  degrees  of  freedom.  Its  similarity  with  the  measured  distributions  suggests  that  the  magnitude  errors  computed  by  the  photometric 
pipeline  are  reliable  (they  may  be  slightly  underestimated,  by  about  10%,  in  the  u  band).  [See  the  electronic  edition  of  the  Journal  for  a  color  version  of  this  figure .] 


The  spatial  variation  of  the  internal  zero-point  error  can  be  sep¬ 
arated  into  “color”  errors,  relative  to  a  fiducial  band,  say  r,  and  an 
overall  “gray”  error  (e.g.,  unrecognized  temporal  changes  in  at¬ 
mospheric  transparency  due  to  gray  clouds), 


Below,  we  discuss  methods  for  estimating  both  the  gray  error 
<5,.(R.A.,  decl.)  and  the  color  errors  <5ffl,-(R.A.,  decl.). 

The  deviation  of  the  internal  SDSS  system  from  a  perfect  AB 
system,  A  m ,  can  also  be  expressed  relative  to  the  fiducial  r  band, 


<5m(R.A.,decl.)=  <Sr(R.A.,  decl.)+  <W(R-A.,  decl.).  (5) 


A  m  —  Ar+  Amr. 


(6) 
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Fig.  2. — Comparison  of  the  color-magnitude  and  color-color  distributions  for  variable  and  nonvariable  unresolved  sources  from  SDSS  stripe  82.  The  distributions  of 
sources  with  an  rms  scatter  of  the  g- band  magnitude  below  0.05  mag  are  shown  by  linearly  spaced  contours  (nonvariable  sources).  Sources  with  320°  <  R.  A.  <  330°,  with 
rms  scatter  in  the  g  band  larger  than  0.05  mag  and  y2  greater  than  3,  are  shown  by  dots  (variable  sources).  The  dots  are  color-coded  according  to  the  observed  rms  scatter  in 
the  g  band  (0.05-0.10  mag;  see  the  legend,  where  red  indicates  variability  larger  than  0.1  mag).  Note  how  low-redshift  quasars  (u  —  g  <  0.6)  and  RR  Lyrae  stars 
(u  —  g  ~  1.1,  g  —  r  ~  0;  Ivezic  et  al.  2005)  clearly  stand  out  as  variable  sources  ( red  dots). 


The  motivation  for  this  separation  is  twofold.  First,  Amr  can  be 
constrained  by  considering  the  colors  (spectral  energy  distri¬ 
butions)  of  objects  independently  from  the  overall  flux  scale  (this 
can  be  done  using  both  external  observations  and  models).  Sec¬ 
ond,  it  is  difficult  to  find  a  science  result  that  crucially  depends  on 
knowing  the  “gray-scale”  offset,  Ar,  at  the  l%-2%  level.  On  the 
other  hand,  knowing  the  “band-to-band”  offsets,  Amr,  with  such 
an  accuracy  is  important  for  many  applications  (e.g.,  photometric 
redshifts  of  galaxies,  Type  la  supernova  cosmology,  and  testing  of 
stellar  and  galaxy  models). 

Fitting  the  SDSS  spectra  of  hot  white  dwarfs  to  models, 
Eisenstein  et  al.  (2006 )  determined  the  AB  color  corrections  A  mr 
to  be  —0.040,  +0.000,  +0.015,  and  +0.030  mag  for  m  =  ugiz, 
respectively,  with  an  uncertainty  of  ^0.01-0.02  mag.  Similar 


values  of  Amr,  also  based  on  hot  white  dwarfs,  were  obtained 
by  Holberg  &  Bergeron  (2006):  -0.046,  +0.001,  +0.013,  and 
+0.024  mag  for  m  =  ugiz,  respectively.  It  may  be  possible  to 
determine  these  corrections  with  an  uncertainty  of  ^0.01  mag, 
and  such  efforts  are  in  progress  (J.  Marriner  2007,  private 
communication). 

The  overall  gray  flux  scale  calibration  error,  Ar,  is  determined 
by  the  accuracy  of  the  absolute  flux  calibration  of  fundamental 
standard  BD  +17  4708  (Fukugita  et  al.  1996),  the  accuracy  of 
tying  the  primary  standard  star  network  to  BD  +17  4708,  the 
accuracy  of  transferring  the  primary  standard  star  network  to  the 
secondary  standard  star  network,  and  the  accuracy  of  the  cali¬ 
bration  of  the  survey  imaging  data  using  the  secondary  standard 
star  network.  Given  these  numerous  sources  of  error,  it  seems 
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Fig.  3 . — Color-magnitude  and  color-color  distributions  of  candidate  standard  stars  brighter  than  g  —  20.  The  dots  are  color-coded  according  to  the  observed  rms  scatter 
in  the  g  band  (0-0.05  mag;  see  the  legend).  Note  the  absence  of  quasars  and  RR  Lyrae  stars  visible  in  Fig.  2. 


unlikely  that  Ar  <  0.02  mag.  On  the  other  hand,  formal  anal¬ 
ysis  of  all  the  error  contributions,  as  well  as  a  direct  comparison 
to  Hubble  Space  Telescope  (1 1ST  )  observations  of  hot  white 
dwarfs  (Eisenstein  et  al.  2006),  suggests  that  Ar  does  not  ex¬ 
ceed  0.05  mag;  in  addition,  the  analysis  of  HST spectrophotometry 
and  other  ground-based  data  by  Holberg  &  Bergeron  (2006)  yields 
A  r  =  0 . 003  mag.  N  ote,  however,  that  all  these  uncertainties  in  the 
definition  and  transfer  of  the  standard  star  network  become  moot 
if  one  accepts  the  following: 

1 .  We  do  not  need  to  know  A exquisitely  well  for  most  scien¬ 
tific  applications.30  Even  if  we  do,  this  is  just  a  single  number  that 
modifies  the  cataloged  photometry  for  all  the  sources  and  all  the 
bands  in  the  same  fashion. 

2.  Uncertainties  in  the  determination  of  Amr  are  of  the  order 
of  0.01  mag. 

3.  We  can  constrain  or  correct  for  £m(R.A.,  deck)  at  the 
0.01  mag  level. 


30  Here  we  assume  that  color  offsets  between  different  surveys  and  photomet- 
ric  systems  can  be  treated  in  a  fashion  similar  to  that  of  the  color  offsets  between 
SDSS  bands.  For  example,  hot  white  dwarfs  could  be  used  to  compute  the  AB 
color  correction,  Amr,  between  the  SDSS  z  band  and  2MASS  J  band  with  an  ac¬ 
curacy  of  ~0. 01-0.02  mag. 


In  other  words,  the  band-to-band  calibration  can  be  fixed  by  adopt¬ 
ing  Amr,  determining  and  correcting  for  <5,„(R.A.,  deck)  guaran¬ 
tees  internal  consistency,  and  the  only  remaining  relatively  free 
parameter  is  A,..  Such  a  system  is  then  no  longer  defined  by  a  set 
of  celestial  standards  but  rather  by  the  functions  fb-  The  catalog 
presented  here  is  one  realization  of  such  a  photometric  system. 

2.5.1.  Determination  of  <5,„(R.A.,  deck) 

We  now  proceed  to  describe  methods  for  constraining  t)m(  R.A., 
deck).  The  region  covered  by  SDSS  stripe  82  is  an  elongated 
rectangle  with  an  aspect  ratio  of  1 : 50  and  with  the  long  side  par¬ 
allel  to  the  celestial  equator.  Because  of  this  large  aspect  ratio, 
and  because  different  effects  contribute  to  the  right  ascension  and 
declination  dependences  of  6m,  we  assume  that  it  can  be  expressed 
as  a  sum  of  independent  functions  of  right  ascension  or  declina¬ 
tion,  respectively, 

UR.A,  deck)  =  «f(decl.)+  6?(R.A.),  (7) 

with  (<5®(decl.)}RA  =  0  and  (<5“l(R.A.))dcc|  =  0,  where  (S)x 
denotes  the  average  of  S  over  direction  x. 

The  first  term,  <5^(deck),  is  dominated  by  the  errors  in  the 
flat-field  vectors  (for  drift  scanning,  flat-field  corrections  are 
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Fig.  4. — Same  as  Fig.  3,  but  only  for  variable  sources  and  with  a  different  color  coding  (0.05-0.10  mag  range  of  the  g-band  rms,  instead  of  0-0.05  mag). 


one-dimensional).  The  flat-field  determination  for  SDSS  was 
difficult  due  to  scattered  light,31  and  there  are  probably  systematic 
errors  in  the  stellar  photometry  at  the  0.01  mag  level  in  the  griz 
bands  and  the  0.02  mag  level  in  the  u  band  (perhaps  somewhat 
larger  at  the  edge  of  the  imaging  camera),  as  we  demonstrate  be¬ 
low.  Since  these  systematic  errors  do  not  cancel  when  averaging 
many  observing  runs  (because  most  stars  are  always  observed  by 
the  same  CCD  and  fall  on  roughly  the  same  position  within  the 
CCD),  <5®(decl.)  could  be  as  large  as  ~0. 01-0. 02  mag  on  spatial 
scales  much  smaller  than  the  chip  width  of  13'. 

The  second  term,  7“l(  R.A.),  is  dominated  by  unrecognized 
fast  variations  of  atmospheric  extinction  (e.g.,  due  to  cirrus),  be¬ 
cause  for  each  observing  run  only  a  single  zero  point  per  CCD  is 
determined.32  While  such  variations  are  uncorrelated  for  differ¬ 
ent  runs,  it  is  possible  that  they  do  not  average  out  fully  at  the  1  % 
level. 


31  For  details,  see  http://www.sdss.org/dr5/algorithms/flatfield.html  and 
Stoughton  et  al.  (2002). 

32  Although  not  relevant  for  the  equatorial  runs  discussed  here,  the  slow 
change  of  atmospheric  extinction  due  to  varying  air  mass  is  accounted  for;  the 
response  of  the  telescope  and  CCDs  is  stable  at  the  <1%  level  on  single-night 
timescales. 


2.5.2.  The  Six  ugriz  Photometric  Systems 

Before  we  describe  methods  for  determining  <5m(R.A.,  decl.), 
we  address  the  bandpass  differences  between  the  six  camera  col¬ 
umns.  The  bandpasses  had  been  measured  for  each  CCD  using  a 
monochromator  (M.  Doi  et  al.  2007,  in  preparation)  and  found 
not  to  be  identical,  as  shown  in  Figure  6.  These  differences  be¬ 
tween  the  bandpasses  induce  color-tenn  errors  in  the  reported 
SDSS  photometry  because  the  magnitudes  of  calibration  stars 
obtained  by  the  PT  are  transformed  to  the  2.5  m  system  using  a 
single  set  of  color  terms.33  In  other  words,  the  color  difference 
between  a  blue  star  and  a  red  star  depends  on  which  camera  col¬ 
umn  the  stars  fall  on;  at  the  ~7  %  accuracy  level  discussed  here, 
there  are  six  SDSS  ugriz  systems. 

We  used  the  measured  response  curves  to  generate  synthetic 
ugriz  photometry  corresponding  to  six  camera  columns  (via  eq.  [1]) 
for  175  stars  from  the  atlas  of  Gunn  &  Stryker  (1983).  The  dif¬ 
ferences  between  the  predicted  magnitudes  for  each  camera  col¬ 
umn  and  the  values  generated  with  the  response  curves  which 
define  the  SDSS  system34  represent  photometric  corrections  due 


33  See  http://www.sdss.org/dr5/algorithms/jeg_photometric_eq_drl .html. 

34  Transmission  curves  for  the  SDSS  2.5  m  photometric  system  are  available 
at  http://www.sdss.org/dr5/instmments/imager. 
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Fig.  5. — Top :  Sky  density  of  all  the  point  sources  with  at  least  four  observa¬ 
tions  in  each  of  the  g,  r,  and  i  bands  (solid  line)  and  of  selected  candidate  standard 
(nonvariable)  stars  (dashed  line).  Bottom :  Ratio  of  the  two  curves  shown  in  the 
top  panel.  For  reference,  Galactic  coordinates  (l,  b)  are  (46,  —24),  (96,  —60),  and 
(190,  —37)  for  o;j2ooo.o  —  —50°,  0°,  and  60°  (at  <5j2000.0  =  0°).  [See  the  electronic 
edition  of  the  Journal  for  a  color  version  of  this  figure .] 

to  bandpass  differences.  These  corrections  may  be  as  large  as 
0.02  mag  for  the  reddest  stars  but  fortunately  admit  simple  lin¬ 
ear  fits  as  a  function  of  color  (u  —  g  for  the  u  band,  g  —  r  for  g 
and  r,  and  r  —  i  for  i  and  z).  The  u  band  is  the  only  one  where 
piecewise  fits  are  required  (in  the  range  0.7  <u  —  g  <  1.3,  the 
g  —  r  color  provides  a  better  fit  than  the  u  —  g  color).  We  have 
applied  these  best-fit  corrections  to  all  sources  in  the  catalog. 
The  median  and  rms  scatter  for  the  distributions  of  corrections 
evaluated  for  all  stars,  and  for  each  color  and  camera  column 
combination,  are  listed  in  Table  1 . 

The  bandpass  differences  have  the  largest  impact  on  the  i  —  z 
color  of  red  sources.  The  color-term  errors  result  in  a  rotation  of 
the  stellar  locus  in  the  i  —  z  versus  r  —  i  color-color  diagram  (see 
Fig.  2,  bottom  right),  and  we  use  this  fact  to  demonstrate  the  im¬ 
provement  in  photometry  due  to  applied  corrections.  We  use  the 
mean  position  of  the  stellar  locus  to  “predict”  the  i  —  z  color  from 
the  measured  r  —  i  color  and  compute  the  difference  between  pre¬ 
dicted  and  measured  i  —  z  colors  separately  for  blue  (0.1  <  r 
i  <  0.2)  and  red  (0.8  <  r  —  i  <  1.4)  stars  and  in  small  bins  of 
declination  (cross-scan  direction).  The  difference  of  these  resid¬ 
uals  (blue  vs.  red  stars)  effectively  measures  the  locus  position 
angle  and  is  not  sensitive  to  photometric  zero-point  errors  and 
flat-field  errors  (which  can  only  induce  locus  shifts  that  have  no 
effect  on  this  test  because  they  cancel  out). 

The  top  panel  in  Figure  7  shows  the  median  i  —  z  residuals  be¬ 
fore  applying  corrections  for  different  transmission  curves,  and 
the  bottom  panel  shows  the  results  based  on  corrected  photometry. 
The  rms  scatter  decreases  from  9  to  3  mmag  after  applying  the 
corrections.  The  remaining  deviations  of  residuals  from  zero  could 
be  due  to  the  fact  that  we  have  not  corrected  for  the  dependence 
of  the  flat  fields  on  source  color  (measurements  of  this  depen¬ 
dence  are  not  available).  Because  this  is  currently  the  only  avail¬ 
able  SDSS  catalog  with  photometry  corrected  for  bandpass 
differences,  it  effectively  represents  a  practical  definition  of  the 
SDSS  photometric  system  (i.e.,  it  provides  photometry  for  real 
sources  on  the  sky;  the  system  is  formally  defined  by  6b,  of 
course). 


Fig.  6. — Top :  Comparison  of  the  five  transmission  curves,  0  (see  eq.  [3]  for 
definition),  that  define  the  SDSS  photometric  system  (dashed  lines)  to  measured 
transmission  curves  sorted  by  camera  column  and  shown  by  solid  lines,  as  marked 
in  the  panel  (the  area  under  each  curve  is  unity  by  definition).  Bottom :  Differences 
between  the  measured  transmission  curves  and  the  curves  that  define  the  system. 
These  differences  induce  color  terms  that  result  in  systematic  photometric  errors 
as  a  function  of  source  color.  The  largest  color  terms  are  present  for  the  z  band  in 
camera  columns  2  and  6,  with  the  errors  well  described  by  Az2  =  Z2  —  zsdss  — 
— 0.019(r  —  i)  and  Az6  =  +0.01  l(r  —  i)  mag,  respectively. 

2.5.3.  Other  Sources  of  Systematic  Errors 

The  nonlinearity  of  the  detectors  (as  a  function  of  source  bright¬ 
ness)  has  also  been  measured  in  situ  by  M.  Doi  et  al.  and  found  to 
be  at  most  a  couple  of  percent  effect  over  the  relevant  dynamic 
range.35  These  corrections  are  determined  with  a  sufficient  accu¬ 
racy  (<5  mmag  impact  on  photometry)  and  are  already  imple¬ 
mented  in  the  photometric  pipeline  (i.e.,  before  performing 
photometric  calibration). 

Similarly  to  the  camera  column-to-column  bandpass  differ¬ 
ences,  variations  in  the  wavelength  dependence  of  atmospheric 
transmissivity  can  also  induce  systematic  errors  that  depend  on 
source  colors.  Assuming  a  standard  atmosphere,  using  synthetic 
photometry  for  stars  from  the  Gunn-Stryker  atlas  we  find  that 
this  effect  can  induce  offsets  of  up  to  ~0.01  mag  for  the  u  —  g  and 
g  —  r  colors  when  air  mass  is  varied  by  0.3  from  its  fiducial  value 
of  1.3.  Flowever,  because  all  stripe  82  data  are  obtained  at  the 
same  air  mass,  this  effect  is  not  relevant  for  the  catalog  discussed 
here.  At  least  in  principle,  similar  and  potentially  larger  errors 
could  be  induced  even  at  a  constant  air  mass  if  the  wavelength 
dependence  of  atmospheric  transmissivity  is  significantly  differ¬ 
ent  from  the  assumed  standard  atmosphere.  Given  that  such  er¬ 
rors  would  probably  average  out,  and  that  there  are  no  available 
measurements  of  the  wavelength  dependence  of  atmospheric 
transmissivity36  for  the  SDSS  data  considered  here,  we  ignore 
this  effect  hereafter. 

2.6.  Determination  of  Flat-Field  Corrections 

We  use  two  methods  based  on  SDSS  data  to  constrain 
h,„(R.A.,  deck):  a  direct  comparison  with  the  secondary  standard 


35  For  details,  see  http://www.sdss.org/dr5/instraments/imager/nonlinearity 
.html. 

36  Some  handle  on  the  stability  of  the  wavelength  dependence  of  atmospheric 
transmissivity  can  be  obtained  by  studying  first-order  extinction  coefficients  de¬ 
termined  by  the  photometric  calibration  pipeline.  They  show  a  cyclic  variation  dur¬ 
ing  the  year,  with  an  nns  scatter  of  residuals  around  the  mean  relation  of  0.01  mag. 
The  amplitude  of  the  yearly  variation  of  the  /  -band  extinction  coefficient  is  —20% 
about  the  mean  value,  and  the  wavelength  dependence  of  the  variation  appears 
consistent  with  the  addition  of  a  gray  opacity  source  during  the  summer  months. 
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TABLE  1 

The  Color-Term  Corrections 


Column  (du)  adu  ( dg )  adg  {dr)  adr  ( di )  adi  {dz)  ad- 


1  .  1.0  0.3  1.8  0.4  -1.2  0.3  0.0  0.2  0.6  0.5 

2  .  11.0  2.3  7.6  1.6  -1.1  0.2  -0.0  0.2  -0.8  1.3 

3  .  -5.4  0.8  -2.0  0.0  2.9  1.0  -0.2  0.3  -0.2  0.5 

4  .  -9.5  2.1  -1.0  0.2  -0.9  0.1  0.0  0.0  -0.4  0.3 

5  .  3.9  1.0  -4.4  1.3  -0.0  0.2  -0.2  0.1  0.4  0.7 

6  .  -0.5  0.8  -2.6  0.4  -0.6  0.1  -0.2  0.3  0.9  1.2 


Note. — The  median  and  rms  scatter  for  photometric  corrections  that  place  the  measurements  in  six  camera  columns  on  the  survey  sys¬ 
tem  (in  millimagnitudes). 


star  network  (§  2.6.2)  and  a  method  based  on  stellar  colors  (§2.6.1). 
Each  of  these  methods  has  its  advantages  and  disadvantages.  We 
perform  tests  of  the  final  catalog  quality  using  a  method  based  on 
the  photometric  redshift  relation  for  galaxies  (§  2.7.2),  and  we 
also  compare  the  SDSS  photometry  to  an  independent  set  of  stan¬ 
dards  provided  by  Stetson  (2000,  2005)  in  §  2.7.3. 

A  determination  of  the  hm  ( R .  A . .  deck)  error  from  a  direct  com¬ 
parison  with  the  secondary  standard  star  network  (hereafter  “PT 
comparison”)  might  be  considered  the  best  method  a  priori.  How¬ 
ever,  it  is  quite  possible  that  the  secondary  standard  star  network 
itself,  from  which  this  catalog  is  calibrated,  may  induce  a  spatial 
variation  of  the  photometric  zero  points  at  the  0.01  mag  level 
(Smith  et  al.  2002).  In  addition,  there  are  not  enough  stars  to  con¬ 
strain  <5;^(decl.)  with  sufficient  spatial  resolution  (say,  at  least 
~100  pixels,  or  ^0.01°).  For  example,  there  are  ^20,000  sec¬ 
ondary  standards  from  stripe  82  in  the  averaged  catalog  that  are 
not  saturated  in  the  gri  bands  in  the  2.5  m  scans  and  have  PT  er¬ 
rors  smaller  than  0.03  mag  (^8000  stars  are  usable  in  the  z  band 
and  only  ^3000  in  the  u  band).  If  these  stars  are  binned  in  the  dec¬ 
lination  direction  every  0.01°  (250  bins,  each  90  pixels  wide),  b® 
in  each  bin  can  be  constrained  to  about  ^0.005  mag  (0.01  mag  in 
the  u  band).  This  is  barely  sufficient  in  the  gri  bands  and  cannot 


Dec  (deg) 

Fig.  7. — Test  based  on  the  stellar  locus  position  angle  in  the  i  —  z  vs.  r  —  i  color- 
color  diagrams  that  demonstrates  the  existence  of  color  terms  between  different 
camera  columns.  The  y-axis  shows  the  difference  in  the  i  —  z  color  residuals  of 
blue  (0.1  <  r  —  i  <  0.2)  and  red  (0.8  <  r  —  i  <  1 .4)  stars,  with  the  residuals  com¬ 
puted  as  the  difference  between  measured  i  —  z  colors  and  those  predicted  using 
the  mean  stellar  locus  (see  Fig.  3,  bottom  right  panel).  The  vertical  solid  lines  mark 
the  approximate  boundaries  between  different  camera  columns,  with  the  vertical 
dashed  lines  marking  approximate  boundaries  between  the  “north”  and  “south” 
strips  in  a  stripe  (see  §2.1  for  definitions).  The  top  panel  shows  the  results  before 
applying  corrections  for  different  transmission  curves  (see  Fig.  6),  and  the  bottom 
panel  shows  the  results  based  on  corrected  photometry.  It  is  evident  that  the  re¬ 
siduals  in  the  bottom  panel  are  much  smaller,  with  rms  scatter  decreasing  from  9  to 
3  mmag.  [See  the  electronic  edition  of the  Journal for  a  color  version  of  this figure.] 


provide  satisfactory  constraints  on  the  flat-fielding  errors  in  the  u 
band  (where,  unfortunately,  these  errors  are  the  largest).  Similarly, 
^“‘(R.A.)  can  be  constrained  in  0.5°  wide  right  ascension  bins 
with  a  similar  accuracy,  but  the  secondary  standard  stars  are  not 
uniformly  distributed  in  right  ascension.  For  these  reasons,  we 
combine  the  PT  comparison  with  the  stellar  locus  method  to  de¬ 
termine  flat-field  corrections. 

2.6.1.  Color  Corrections  from  the  Stellar  Locus  Method 

The  stellar  distribution  in  color-color  space  at  high  Galactic  lat¬ 
itudes37  (\b\  >  30)  is  remarkably  uniform  at  the  faint  flux  levels 
probed  by  SDSS,  as  discussed  in  detail  by  Ivezic  et  al.  (2004a). 
Systematic  photometric  errors,  other  than  an  overall  gray  error, 
manifest  themselves  as  shifts  in  the  position  of  the  stellar  locus 
that  can  be  tracked  using  the  four  principal  colors  (sw'.rv)  defined 
by  Ivezic  et  al.  (2004a).  These  colors  are  linear  combinations  of 
magnitudes, 

P2  —  riti  -4-  Bg  -f-  Cr  -p  Di  -P  Ez  -P  F ,  (8) 

where  P2  =  s ,  w,  x,y,  and  measure  the  distance  from  the  center 
of  the  locus  in  various  two-dimensional  projections  of  the  four¬ 
dimensional  stellar  color  distributions  (s,  perpendicular  to  the  blue 
part  of  the  locus  in  the  g  —  r  vs.  u  —  g  plane;  w,  perpendicular  to 
the  blue  part  in  the  r  —  i  vs.  g  —  r  plane;  x,  perpendicular  to  the  red 
part  with  g  —  r  ~  1 .4  in  the  r  —  i  vs.  g  —  r  plane;  and  y,  perpen¬ 
dicular  to  the  locus  in  the  i  —  z  vs.  r  —  i  plane).  The  matrix  of  co¬ 
efficients  A-F  is  listed  in  Table  2  (for  more  details,  see  Ivezic  et  al. 
2004a).  Of  course,  the  measurements  must  be  corrected  for  the 
effects  of  interstellar  dust  extinction;  we  use  maps  provided  by 
Schlegel  et  al.  (1998,  hereafter  SFD98).  The  properties  of  two  of 
these  colors  (s  and  w)  are  illustrated  in  Figure  8. 

The  fact  that  the  median  principal  colors  are  close  to  zero  shows 
that  the  averaging  procedure  did  not  induce  any  shifts  in  zero 
points  compared  to  the  average  of  291  SDSS  runs  which  were 
used  to  define  the  principal  colors  (Ivezic  et  al.  2004a).  The  same 


37  At  low  Galactic  latitudes  several  effects,  discussed  below,  prevent  the  use 
of  this  method  for  calibration  purposes. 

TABLE  2 

The  swxy  Principal  Color  Definitions 


PC  A  B  C  D  E  F 


s .  -0.249  0.794  -0.555  0.0  0.0  0.234 

w .  0.0  -0.227  0.792  -0.567  0.0  0.050 

x .  0.0  0.707  -0.707  0.0  0.0  -0.988 

j .  0.0  0.0  -0.270  0.800  -0.534  0.054 


No.  3,  2007 


SDSS  STANDARD  STAR  CATALOG  FOR  STRIPE  82 


983 


s  w 

Fig.  8. — Top  row :  Two  projections  of  the  stellar  locus  (rotated  locus  from  the  ugr  and  gri  planes;  see  Fig.  3)  constructed  using  averaged  photometry.  The  large  dots 
show  the  median  values  of  the  s  and  w  principal  colors  (perpendicular  to  the  locus  at  its  blue  edge)  in  bins  of  the  principal  color  along  the  locus.  Middle  and  bottom  rows : 
Histograms  for  each  principal  color  on  linear  and  logarithmic  scales  (essentially  the  locus  cross  sections).  The  dark  histograms  are  constructed  using  the  averaged  pho¬ 
tometry,  and  the  light  histograms  are  based  on  single-epoch  photometry.  The  former  are  narrower  than  the  latter  due  to  increased  photometric  accuracy.  The  best-fit  Gaussians, 
with  parameters  listed  in  the  middle  row,  are  shown  by  the  dashed  lines  in  the  bottom  row  (only  for  the  averaged  photometry).  [See  the  electronic  edition  of  the  Journal  for  a 
color  version  of  this  figure .] 


conclusion  is  reached  by  comparing  the  averaged  photometry  with 
the  secondary  standard  star  network:  the  median  photometric  re¬ 
siduals  at  the  so-called  crossing  colors38  are  4, 6, 3, 2,  and  2  mmag 
in  the  ugriz  bands,  respectively.  Yet  another  test  is  a  direct  compar¬ 
ison  of  averaged  photometry  with  single-epoch  photometry.  Using 
the  SDSS  Data  Release  5  photometry,  we  find  that  the  largest  me¬ 
dian  magnitude  difference  between  the  two  sets  is  2  mmag  in  the 
u  band. 

It  is  noteworthy  that  the  widths  of  the  principal  color  distri¬ 
butions  (i.e.,  the  thickness  of  the  stellar  locus)  constructed  with 
averaged  photometry  are  much  smaller  than  when  using  single¬ 
epoch  data  (see  the  bottom  four  panels  in  Fig.  8).  Indeed,  all  four 
principal  color  distributions  are  “resolved”  using  this  high- 
quality  photometry  (see  the  fifth  column  in  Table  3). 


38  Crossing  colors  are  roughly  the  median  colors  of  the  observed  stellar  pop- 
ulation;  for  details,  see  http://www.sdss.org/dr5/algorithms/jeg_photometric_ 
eq_drl.html. 


TABLE  3 

The  Effect  of  Repeated  Measurements  on  the  Width  of  the  Stellar  Locus 


Width  for 

PC  rms  1  Obs.a  rms  N  Obs.b  Median0  PC/PC  Errord 


5 .  31  19  3.0  3.02 

w .  25  10  1.1  1.82 

x .  42  34  1.2  5.34 

y .  23  9  0.8  1.64 


a  The  locus  width  determined  using  single-epoch  SDSS  observations  (in 
millimagnitudes). 

b  The  locus  width  determined  using  multiple  SDSS  observations  (in 
millimagnitudes). 

c  The  median  principal  color  determined  using  multiple  SDSS  observations 
(in  millimagnitudes). 

d  The  locus  width  normalized  by  expected  measurement  errors. 
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TABLE  4 

The  Flat-Field  Corrections 

Vol.  134 

Band 

Width3 

Minimum1* 

Maximum0 

ur . 

.  22 

-53 

53 

gr . 

.  12 

-27 

19 

r . 

.  7 

-17 

17 

ri . 

.  4 

-10 

13 

rz . 

.  7 

-14 

19 

Note. — The  r-band  correction  is  determined  using  observations  with  the  PT, 
and  the  ugiz  corrections  are  determined  using  the  stellar  locus  method  (see  §§  2.6. 1 
and  2.6.2). 

a  The  rms  scatter  for  applied  flat-field  corrections  (in  millimagnitudes). 
b  The  minimum  value  of  the  applied  correction  (in  millimagnitudes). 
c  The  maximum  value  of  the  applied  correction  (in  millimagnitudes). 


Dec  (deg) 

Fig.  9. — Top  panel :  Distribution  of  the  residuals  between  the  PT  photometry 
and  averaged  magnitudes  in  the  gri  bands.  Second  panel :  Applied  flat-field  cor¬ 
rection  in  the  r  band,  which  was  derived  as  the  mean  of  the  residuals  shown  in  the 
top  panel.  Third  and  fourth  panels :  Applied  flat-field  corrections  in  the  other  four 
bands,  expressed  relative  to  the  r  band  ( third  panel ,  u;  fourth  panel,  giz ),  which 
were  derived  using  the  stellar  locus  colors.  The  low-order  statistics  for  these  cor¬ 
rections  are  listed  in  Table  4.  [See  the  electronic  edition  of  the  Journal for  a  color 
version  of  this  figure.'] 

Because  the  intrinsic  widths  of  the  principal  color  distributions 
are  so  small,  principal  colors  can  be  used  to  efficiently  track  local 
calibration  problems  using  a  small  number  of  stars,  allowing  a 
high  spatial  resolution.  That  is,  we  require  that  the  locus  not  move 
in  the  multidimensional  color  space.  In  practice,  the  deviations  of 
the  principal  colors  from  zero  can  be  inverted,  using  an  appropri¬ 
ate  closure  relation  (see  the  next  section),  to  yield  flat-field  correc¬ 
tions  (Ivezic  et  al.  2004a).  With  bins  0.01°  wide  in  the  declination 
direction,  or  1°  wide  in  the  right  ascension  direction,  the  flat-field 
corrections  can  be  determined  with  an  accuracy  of  5  mmag  or 
better. 

2.6.2.  Gray  Corrections  from  the  Comparison 
with  the  Secondary  Standard  Star  Network 

The  main  advantage  of  the  stellar  locus  method  is  that  it  can 
constrain  8mr  with  high  spatial  resolution.  However,  it  is  insensi¬ 
tive  to  gray  errors,  parameterized  by  6r  (e.g.,  an  overall  gradient 
of  photometric  zero  points  in  the  declination  direction  that  is  the 
same  in  all  five  bands  would  have  no  effect  on  stellar  colors).  On 
the  other  hand,  the  PT  comparison  can  constrain  br,  but  it  does  not 
provide  enough  spatial  resolution  to  derive  flat-field  corrections, 
especially  in  the  u  and  z  bands.  Therefore,  we  combine  these  two 
methods  to  derive  flat-field  corrections  b^(decl.). 

The  median  differences  between  the  averaged  2.5  m  photo¬ 
metry  and  PT  photometry  for  secondary  standard  stars  in  the  gri 


bands  are  shown  in  Figure  9  ( top  panel).  The  median  differences 
are  computed  for  0.01°  wide  bins  and  then  smoothed  by  a  trian¬ 
gular  filter  [v,  is  replaced  by  0.25(y,_i  +  2j/;  +  y,-+i  )]•  The  resid¬ 
uals  in  all  three  bands  display  similar  behavior  and  imply  about 
0.02  peak-to-peak  variation  between  the  center  and  edges  on  each 
CCD  (resulting  in  about  6  mmag  rms  contribution  to  the  overall 
errors),  as  well  as  an  overall  0.01-0.02  mag  tilt.  These  systematic 
errors  may  be  due  to  imperfect  flat-field  vectors  used  to  reduce  the 
data,  incorrectly  determined  scattered-light  correction  (the  two  are 
somewhat  coupled  in  the  data  reduction  procedures),  or  problems 
in  the  PT  itself  (such  as  the  PT  flat  field). 

At  face  value,  these  residuals  could  be  used  to  correct  the  aver¬ 
aged  2.5  m  photometry  in  each  band  (gri)  separately.  However, 
doing  so  introduces  noise  in  the  stellar  principal  colors  of  about 
5  mmag  (rms)  and  suggests  that  the  differences  in  photometric 
residuals  between  the  three  bands  are  dominated  by  PT  measure¬ 
ment  noise.  On  the  other  hand,  the  2.5  m  versus  PT  residuals  do 
contain  information  about  “gray”  errors  that  cannot  be  deter¬ 
mined  using  stellar  locus.  Hence,  we  take  the  mean  value  of  the 
2.5  m  versus  PT  residuals  in  the  gri  bands  to  represent  the  6r  flat- 
field  correction  and  apply  it  to  the  averaged  2.5  m  photometry  in 
the  r  band.  The  applied  r-band  correction  is  shown  in  Figure  9 
(second panel)  and  has  an  rms  scatter  of  7  mmag  (for  250  bins), 
with  the  largest  correction  less  than  0.02  mag. 

In  the  second  step,  we  use  the  stellar  locus  to  derive  the  8mr 
corrections  in  each  band  (ugiz).  The  derivation  of  these  correc¬ 
tions  is  essentially  identical  to  the  procedure  described  by  Ivezic 
et  al.  (2004a).  Also,  together  with  a  PT-based  8r  correction,  this 
is  essentially  the  same  method  as  used  to  derive  flat-field  correc¬ 
tions  for  the  whole  SDSS  survey.39  In  particular,  we  used  the  same 
closure  relation  (the  stellar  locus  method  gives  four  equations  for 
five  unknowns)  that  is  based  on  averaged  2.5  m  versus  PT  resid¬ 
uals  in  the  gri  bands.  The  resulting  flat-field  corrections,  8mr,  in  the 
ugiz  bands  are  shown  in  the  third  and  fourth  panels  in  Figure  9  and 
summarized  in  Table  4. 

Due  to  low  stellar  counts  and  the  strongest  scattered  light,  the 
u-band  correction  is  expected  to  have  the  largest  noise  (~5— 
10  mmag),  which  is  consistent  with  the  observed  behavior.  It  is 
thus  likely  that  some  of  the  variation  on  scales  of  ^0.01°  is  not 
real.  On  the  other  hand,  it  could  be  argued  that  systematic  errors 
could  actually  be  much  larger  on  even  smaller  spatial  scales  but 
get  averaged  out  in  90  pixel  wide  bins.  However,  in  addition  to 
not  having  a  reason  to  believe  in  such  high  spatial  frequency  ef¬ 
fects  (e.g.,  the  sky  background  does  not  show  any  evidence  for 
them),  no  additional  scatter,  except  the  expected  statistical  noise, 
is  observed  when  the  bin  size  is  decreased  by  a  factor  of  4. 


39 


For  details,  see  http://www.sdss.org/dr5/algorithms/flatfield.html. 
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Fig.  10. — Dependence  of  the  position  of  the  stellar  locus  in  ugriz  color  space, 
as  parameterized  by  the  median  principal  colors  swxy,  as  a  function  of  right  as¬ 
cension.  Close  to  the  edges,  the  median  colors  deviate  significantly  from  zero.  This 
is  caused  by  intrinsic  changes  in  the  stellar  locus  due  to  stellar  population  variations 
and  overestimated  interstellar  extinction  corrections  for  red  stars,  rather  than  cali¬ 
bration  problems.  [See  the  electronic  edition  of  the  Journal  for  a  color  version  of 
this  figure .] 


Last  but  not  least,  it  is  important  to  emphasize  that  these  cor¬ 
rections  are  not  setting  photometric  zero  points,  but  only  correct¬ 
ing  for  variations  in  response  across  each  CCD.  As  discussed 
above,  the  AB  photometric  zero  points,  relative  to  the  fiducial 
r  band,  are  effectively  set  by  adopting  values  for  A,,,,.. 
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Fig.  1 1 . — Top  panel :  Implied  photometric  zero-point  errors  based  on  the  stel¬ 
lar  locus  method  (ugiz  from  the  bottom  to  the  top  at  either  edge).  While  the  im¬ 
plied  errors  are  small  for  —25°  <  R.  A.  <  40°,  they  become  exceedingly  large  out¬ 
side  this  range.  This  is  due  to  problems  with  the  stellar  locus  method  rather  than 
problems  with  calibration.  Second  panel :  Median  r-band  extinction  derived  from 
the  SFD98  maps.  Third  and  fourth  panels :  Median  residuals  between  the  PT  pho¬ 
tometry  and  averaged  magnitudes  in  the  uz  ( third  panel )  and  gri  bands  ( fourth 
panel).  The  low-order  statistics  for  these  residuals  are  listed  in  Table  5.  [See  the 
electronic  edition  of  the  Journal  for  a  color  version  of  this  figure .] 


2.7.  Tests  of  Catalog  Quality 
2.1  A.  The  Internal  Tests 

At  least  in  principle,  the  same  methods  used  to  derive  <5, f  (deck) 
could  be  used  to  derive  R.A.).  However,  in  practice  this  is  not 

possible  for  at  least  two  reasons:  first,  the  right  ascension  distri¬ 
bution  of  secondary  standard  stars  is  not  as  unifomi  as  their  decli¬ 
nation  distribution,  and  second,  the  assumption  of  the  constancy 
of  the  stellar  locus  in  color  space  is  invalid  along  the  “long”  scan 
direction  (discussed  below).  For  these  reasons,  we  only  use  the 
PT  comparison  and  stellar  locus  methods  to  estimate  the  level  of 
internal  zero-point  variations  with  right  ascension,  and  we  do  not 
correct  the  data.  In  the  next  section  we  also  use  another  method, 
based  on  galaxy  colors,  as  an  independent  test  of  catalog  integrity. 

Figure  10  shows  the  median  principal  colors  in  bins  of  right 
ascension.  It  is  evident  that  the  principal  colors  are  close  to  zero 
for  right  ascension  in  the  range  —25°  <  R.  A.  <  40°  but  outside 
this  range  deviate  significantly  from  zero.  This  does  not  neces¬ 
sarily  indicate  problems  with  photometric  calibration,  because  the 
stellar  locus  method  is  expected  to  fail  at  low  Galactic  latitudes  for 


several  reasons.  First,  the  mean  metallicity  of  stars  increases  at 
low  Galactic  latitudes,  and  this  change  may  affect  the  s  and  w  col¬ 
ors.  Second,  at  low  latitudes  red  dwarfs  are  not  behind  the  entire 
dust  screen  measured  by  the  SFD98  maps  (see  Juric  et  al.  [2005] 
for  a  discussion  of  this  point),  and  thus  the  x  color  will  be  biased 
blue  (i.e.,  the  colors  of  red  dwarfs  are  overcorrected  for  the  ISM 
reddening).  And  third,  at  low  latitudes  the  dust  column  increases 
quickly  (see  Fig.  11 ),  and  even  small  errors  in  the  assumed  wave¬ 
length  dependence  of  the  dust  extinction,  or  the  extinction  itself 
as  given  by  the  SFD98  maps,  will  have  noticeable  effects  on  the 
principal  colors.  For  these  reasons  it  seems  plausible  that  the  de¬ 
viations  seen  in  Figure  10  are  not  dominated  by  zero-point  errors. 

This  conclusion  is  supported  by  direct  comparison  of  the  av¬ 
eraged  and  PT  photometry  (Fig.  11).  For  example,  the  largest  me¬ 
dian  photometric  residual  between  the  averaged  catalog  and  PT 
observations  in  the  u  band  is  ~0.02  mag  (see  Table  5),  which  is 
much  smaller  than  the  0. 1  mag  discrepancy  implied  by  the  stellar 
locus  method. 

Table  5  shows  that  the  rms  scatter  of  the  median  photometric 
residuals  (evaluated  in  2°  wide  bins  in  the  right  ascension  direction) 
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TABLE  6 

Residuals  from  the  Mean  Color-Redshift  Relation  for  Red  Galaxies 


Band  Residuals3  Widthb  Minimum0  Maximum11  N°  Residualsf 


u .  -2  9  -18  23  175  5 

g .  6  7  -4  17  647  5 

r .  3  7  -7  10  627  3 

i .  4  7  -10  17  621  2 

z .  1  8  -16  15  286  -2 


a  The  median  value  for  the  bin  medians  (in  millimagnitudes).  There  are  24  bins, 
distributed  inhomogeneously  in  the  right  ascension  direction. 
b  The  rms  scatter  for  the  bin  medians  (in  millimagnitudes). 

0  The  minimum  value  for  the  median  residuals  (in  millimagnitudes). 
d  The  maximum  value  for  the  median  residuals  (in  millimagnitudes). 
e  The  median  number  of  stars  per  bin. 

f  The  median  value  of  the  residuals  for  stars  with  colors  within  0.02  mag  of 
the  crossing  colors  (in  millimagnitudes). 

between  the  averaged  catalog  and  PT  observations  is  <0.01  mag 
in  all  five  bands.  Some  of  that  scatter  must  come  from  the  PT  data 
themselves,  and  thus  the  true  scatter  of  the  photometric  zero 
points  in  the  averaged  catalog  is  probably  even  smaller  than  that 
listed  in  Table  5.  In  addition.  Table  5  shows  that  the  averaged  cat¬ 
alog  and  PT  measurements  are  on  the  same  system  to  within  a  few 
millimagnitudes  (using  the  recommended  photometric  transfor¬ 
mations  between  the  two  telescopes  listed  at  the  SDSS  Web  site; 
see  §  2.2). 

We  have  also  compared  the  catalog  presented  here  to  the  pho¬ 
tometric  reductions  described  by  Padmanabhan  et  al.  (2007).  As 
discussed  in  §  1,  they  determined  flat  fields  by  minimizing  the 
errors  in  relative  photometry  of  multiply  observed  stars  over  the 
whole  survey  region.  Hence,  this  comparison  is  an  essentially 
independent  test  of  the  flat-field  corrections  derived  here,  despite 
the  fact  that  both  catalogs  are  based  on  the  same  observations. 
We  bin  the  photometric  differences  between  the  catalogs  in  0.0 1° 
wide  declination  bins  and  compute  the  median  residual  for  each 
bin  and  band.  The  median  values  of  these  medians  represent 
zero-point  offsets  in  each  band  and  are  equal  to  —7,  2,  —1,  2, 
and  5  mmag  in  the  ugriz  bands,  respectively.  The  rms  scatter  of 
the  median  residuals  reflects  systematic  errors  due  to  flat-field 
errors,  and  we  measure  27,  6,  5,  6,  and  8  mmag  in  the  ugriz 
bands,  respectively.  Except  in  the  u  band,  these  values  indicate 
that  the  systematic  flat-field  errors  are  very  small.  In  the  u  band, 
Padmanabhan  et  al.  (2007)  expect  errors  of  about  0.01  mag, 
and  the  distribution  width  of  the  s  color  implies  about  0.01  mag 
for  the  catalog  discussed  here,  predicting  about  14  mmag  instead 
of  the  measured  27  mmag.  It  is  plausible  that  the  //-band  photom¬ 
etry  may  contain  systematic  errors  unrecognized  by  any  of  the 
methods  discussed  here. 

2.7.2.  The  Tests  of  Catalog  Quality  Based  on  Galaxy  Colors 

The  color  distribution  of  galaxies  is  bimodal  (Strateva  et  al. 
2001;  Yasuda  et  al.  2001;  Baldry  et  al.  2004).  Red  galaxies  have 
an  especially  tight  color-redshift  relation  ( Eisenstein  et  al.  200 1 ), 
with  an  rms  of  0.12  mag  for  the//  —  g  color,  0.05  mag  for  the  g  —  r 
color,  and  0.03  mag  for  the  r  —  i  and  i  —  z  colors  (using  model 
magnitudes).  Deviations  from  the  mean  relations  can  thus  be  used 
to  track  local  calibration  problems.  Of  course,  since  this  is  a  color- 
based  method,  it  can  only  constrain  8mr,  and,  because  red  galaxies 
are  faint  in  the  u  band,  cannot  achieve  high  spatial  resolution  in 
this  band.  Nevertheless,  it  is  a  useful  addition  to  the  stellar  locus 
method  because  it  is  independent  of  the  Milky  Way  structure  and 
secondary  star  network  (although  it  is  sensitive  to  errors  in  the 
ISM  dust  extinction  correction). 


Color  rms  (Decl.)a  \  (Decl.)  rms  (R.A.)b  \  (R  A.) 


u-g .  21  1.5  18  1.3 

g-r .  4  1.3  12  2.6 

.  3  1.3  6  3.5 

i-z .  9  2.9  6  2.9 


Note. — The  rms  widths  of  color-residual  distributions  (in  millimagnitudes) 
and  the  widths  of  distributions  of  residuals  normalized  by  statistical  noise  (y ), 
using  mean  color-redshift  relations  (see  text). 

a  The  rms  for  the  declination  direction,  using  0.025°  wide  bins  (0.1°  for  the 
u  —  g  color). 

b  The  rms  for  the  right  ascension  direction,  using  2°  wide  bins  (5°  for  the  u  —  g 
color). 

We  select  19,377  red  galaxies  with  SDSS  spectra  from  the 
redshift  range  0.02-0.36  using  an  empirical  condition, 

0  <  (g  -  r)  —  0.6  -  2.75  x  redshift  <  0.3,  (9) 

and  determine  their  median  colors  as  a  function  of  redshift  using 
0.01  wide  redshift  bins.  The  residuals  from  the  median  color- 
redshift  relation  are  then  binned  by  declination  to  constrain  6%  and 
by  right  ascension  to  constrain  <Yxt.  The  rms  for  color  residuals 
and  the  widths  of  distributions  of  residuals  normalized  by  sta¬ 
tistical  noise  (based  on  quoted  photometric  errors)  are  listed  in 
Table  6. 

The  residuals  binned  in  the  declination  direction  are  generally 
small  and  consistent  with  statistical  noise.  The  largest  deviations 
from  zero  are  seen  for  the  /'  —  z  color,  with  an  rms  of  9  mmag  and 
maximum  deviation  of  17  mmag  (see  Fig.  12,  top).  Although  the 
rms  is  fairly  small,  it  is  a  factor  of  2.9  larger  than  the  expected 
noise.  The  shape  of  the  i  —  z  residuals  is  similar  to  that  of  the  i  —  z 
residuals  for  stars  discussed  in  §  2.5.2  and  shown  in  Figure  7  (bot¬ 
tom).  The  rms  scatter  for  the  difference  between  the  stellar  and  gal¬ 
axy/  —  z  residuals,  shown  in  Figure  \2(bottom),  is6mmag.  Since 
the  rms  for  the  galaxy  residuals  is  larger  (9  mmag),  it  is  plausible 


Dec  (deg) 


Dec  (deg) 

Fig.  12. — Top'.  Dependence  of  the  median  i  —  z  color  residuals  (with  respect 
to  a  mean  color-redshift  relation)  for  red  galaxies  as  a  function  of  declination.  The 
vertical  solid  lines  mark  the  approximate  boundaries  between  different  camera 
columns,  with  the  vertical  dashed  lines  marking  approximate  boundaries  between 
the  "north”  and  “south”  strips  in  a  stripe.  The  small  rms  scatter  of  only  9  mmag 
demonstrates  that  flat-field  corrections  based  on  the  stellar  locus  position  in  color 
space  are  also  applicable  for  galaxies.  Bottom'.  Difference  between  the  values  shown 
in  the  top  panel  and  the  curve  shown  in  the  bottom  panel  of  Fig.  7.  The  rms  scatter 
for  the  residuals  is  6  mmag.  [See  the  electronic  edition  of  the  Journal  for  a  color 
version  of  this  figure .] 
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Fig.  13. — Comparison  of  systematic  color  errors  implied  by  different  methods. 
Note  that  the  errors  implied  by  the  stellar  locus  method  {jagged  line)  become  very 
large  outside  the  —25°  <  cij2ooo.o  <  40°  range.  As  galaxy  colors  (< triangles )  and  a  di¬ 
rect  comparison  with  the  SDSS  secondary  standard  star  network  {circles)  suggest, 
this  is  due  to  problems  with  the  stellar  locus  method  rather  than  problems  with  cal¬ 
ibration.  [See  the  electronic  edition  of  the  Journal  for  a  color  version  of  this  figure .] 


that  the  same  systematic  effect  dominates  the  remaining  photo¬ 
metric  errors  for  both  stars  and  galaxies.  While  it  is  not  clear  what 
the  cause  of  this  (small)  effect  is,  a  plausible  explanation  is  the 
dependence  of  flat  fields  on  source  color.40 

The  results  for  binning  in  the  right  ascension  direction  are 
shown  in  Figure  13,  where  we  compare  different  methods.  The 
rms  for  implied  color  errors  (with  2°  wide  bins)  from  galaxies  is 
0.006  mag  for  the  r  —  i  and  i  —  z  colors,  0.012  mag  for  the  g  —  r 
color,  and  0.01 8  mag  for  the  u  —  g  color.  The  overall  behavior 
of  red  galaxy  color  residuals  agrees  better  with  the  PT  method 
than  with  color  errors  implied  by  the  stellar  locus  method.  In  par¬ 
ticular,  the  large  errors  outside  the  —25°  <  R.A.  <  40°  range  im¬ 
plied  by  the  latter  method  are  not  consistent  with  red  galaxy  color 
residuals.  On  the  other  hand,  both  red  galaxy  colors  and  the  stellar 
locus  seem  to  show  a  trend  that  the  colors  are  redder  around 
R.A.  =  —10°  than  around  R.A.  =  40°.  The  amplitude  of  this 
effect  varies  from  about  0.02  mag  for  the  g—r  color  to  about 
0.01  mag  for  the  i  —  z  color,  while  the  upper  limit  on  such  a 
slope  implied  by  the  PT  comparison  is  <0.01  mag. 


40  It  is  fair  to  ask  whether  the  applied  flat-field  corrections,  derived  from  stel- 
lar  colors,  are  actually  appropriate  for  galaxies.  They  are,  since  the  i  —  z  color  resid¬ 
uals  for  galaxies  without  any  flat-field  corrections  are  about  twice  as  large  as  those 
shown  in  Fig.  12  (top). 


It  is  not  clear  what  the  cause  of  this  discrepancy  is.  The  obvious 
culprit  is  the  correction  for  interstellar  dust  extinction,  but  the  im¬ 
plied  deviation  is  too  large  to  be  explained  by  any  plausible  errors 
in  the  SFD98  maps.  As  shown  in  Figure  1 1 ,  the  median  extinction 
in  the  r  band  for  the  —  1 0°  <  R.  A.  <  40°  range  is  below  0. 1  mag, 
and  the  resulting  median  correction  for,  e.g.,  the  g  —  r  color  is  be¬ 
low  0.04  mag.  Hence,  to  induce  a  0.02  mag  trend  in  the  g  —  r 
color,  the  SFD98  value  for  the  r-band  extinction,  Ar,  would  have 
to  be  in  error  by  0.05  mag  (the  difference  between  the  values  pro¬ 
vided  for  R.A.  =  —10°  and  R.A.  =  40°).  This  implies  relative 
errors  for  the  SFD98  maps  in  the  range  from  50%  (if  Ar  at  R.A.  = 
—  10°  is  underestimated)  to  100%  (if  Ar  at  R.A.  =  40°  is  overes¬ 
timated),  which  seems  unlikely  (although  not  impossible). 

We  conclude  that  the  PT  comparison  provides  a  good  esti¬ 
mate  of  the  remaining  zero-point  errors  in  the  catalog,  as  listed  in 
Table  5,  but  caution  that  we  do  not  understand  the  above  system¬ 
atic  behavior  of  the  stellar  and  galaxy  colors,  and  that  only  the  PT 
constrains  possible  gray  errors.  In  the  next  section  we  discuss  a 
comparison  to  an  external  data  set  that  supports  this  conclusion. 

2.7.3.  An  External  Test  of  Catalog  Quality 
Based  on  Stetson  s  Standards 

The  only  large  external  data  set  with  sufficient  overlap,  depth, 
and  accuracy  to  test  the  quality  of  the  stripe  82  catalog  is  that  pro¬ 
vided  by  Stetson  (2000, 2005 ).  Stetson’s  catalog  lists  photometry 
in  the  B  VR1  bands  (Stetson’s  photometry  is  tied  to  Landolt’s  stan¬ 
dards)  for  ~1 200  stars  in  common  (most  have  V  <  19.5).  We 
synthesize  the  B VRI  photometry  from  the  SDSS  gri  measure¬ 
ments  using  photometric  transformations  of  the  following  form: 

OTStetson  ~  tn sdss  =  Ac*  +  Bc~  +  Cc  +  D,  (10) 

where  mstetson  =  (BVR1),  otsdss  =(g,g,  r,  i),  and  the  color  c  is 
measured  by  SDSS  ( g  —  r  for  the  B  and  V  transformations  and 
r  —  i  for  the  R  and  I  transformations).  The  measurements  are  not 
corrected  for  ISM  reddening.  Traditionally,  such  transformations 
are  assumed  to  be  linear  in  color.41  We  use  higher  order  terms  in 
equation  (10)  because  at  the  1  %-2%  level  there  are  easily  detect¬ 
able  deviations  from  linearity  for  all  color  choices  (for  details  and 
plots,  see  Ivezic  et  al.  2007).  The  best-fit  coefficients  for  the  trans¬ 
formation  of  SDSS  gri  measurements  to  the  BVR1  system42  and 
low-order  statistics  for  the  mstetson  —  «sdss  difference  distribu¬ 
tion43  are  listed  in  Table  7.  We  find  no  trends  as  a  function  of  mag¬ 
nitude  at  the  <0.005  mag  level.  With  the  listed  transformations, 
the  SDSS  catalog  described  here  could  also  be  used  to  calibrate 
the  data  to  the  B  VRI  system  with  a  negligible  loss  of  accuracy  due 
to  transformations  between  the  two  systems. 

The  B  VRI  photometry  from  Stetson  and  that  synthesized  from 
SDSS  agree  at  the  level  of  0.02  mag  (rms  scatter  for  the  magni¬ 
tude  differences  of  individual  stars;  note  that  the  systems  are  tied 
to  each  other  to  within  a  few  millimagnitudes  by  the  transforma¬ 
tions  listed  in  Table  7).  This  scatter  is  consistent  with  the  claimed 
accuracy  of  both  catalogs  (the  magnitude  differences  normalized 
by  the  implied  error  bars  are  well  described  by  Gaussians  with 
widths  in  the  range  0.7-0. 8).  This  small  scatter  allows  us  to  test 


41  For  various  photometric  transformations  between  the  SDSS  and  other 
systems,  see  Abazajian  et  al.  (2005)  and  http://www.sdss.org/dr5/algorithms/ 
sdssUBVRITransform.html. 

42  The  same  transformations  can  be  readily  used  to  transform  measurements 
in  the  BVR1  system  to  the  corresponding  gri  values  because  B  —  V  =  f(g  —  r) 
and  R  —  I  =  f(r  —  i)  are  mono  tonic  functions. 

43  Note  that  these  transformations  are  valid  only  for  main-sequence  stars  with 
colors  in  the  range  g  —  r  >  0.2  and  r  —  i  <1.5  (roughly,  0.3  <  B  —  V  <  1.6). 
Extrapolation  outside  this  range  may  result  in  large  errors. 
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TABLE  7 

Comparison  with  Stetson’s  Standards:  gri-TO-BVRI  Transformations 


Color  Residuals  mc(ja  crmedb  Xmed°  Residualsalld  craUe  Af  B f  Cf  Df 


B-g .  -1.6  8.7  1.4  1.0  32  0.2628  -0.7952  1.0544  0.0268 

V-g .  0.8  3.9  1.0  0.9  18  0.0688  -0.2056  -0.3838  -0.0534 

R-r .  -0.1  5.8  0.9  1.2  15  -0.0107  0.0050  -0.2689  -0.1540 

I-i .  0.9  6.1  1.0  1.2  19  -0.0307  0.1163  -0.3341  -0.3584 


Notes. — These  transformations  are  valid  only  for  main-sequence  stars  with  colors  in  the  range  g  —  r  >  0.2  and  r  —  i  <  1.5  (roughly,  0.3  <  B  —  V  <  1.6).  Extrapo¬ 
lation  outside  this  range  may  result  in  large  (0.1  mag  or  larger  for  hot  white  dwarfs)  errors. 

a  The  median  value  of  median  transformation  residuals  (differences  between  the  measured  values  of  colors  listed  in  the  first  column  and  those  synthesized  using  eq.  [  1 0])  in 
0.1  mag  wide  g  —  r  bins  for  stars  with  0.25  <  g  —  r  <  1.45  (in  millimagnitudes).  These  medians  of  medians  measure  the  typical  level  of  systematics  in  the  gri-to-B  VRI  pho¬ 
tometric  transformations  introduced  by  the  adopted  analytic  form  (see  eq.  [10]). 
b  The  rms  scatter  for  median  residuals  described  above  (in  millimagnitudes). 

c  The  rms  scatter  for  residuals  normalized  by  statistical  noise.  The  listed  values  are  ~  1 ,  which  indicates  that  the  scatter  around  adopted  photometric  transformations 
discussed  in  footnote  a  is  consistent  with  expected  noise. 

d  The  median  value  of  residuals  evaluated  for  all  stars  (in  millimagnitudes). 

e  The  rms  scatter  for  residuals  evaluated  for  all  stars  (in  millimagnitudes). 

f  Coefficients  A—D  needed  to  transform  SDSS  photometry  to  the  BVR1  system  (see  eq.  [10]). 


for  the  spatial  variation  of  zero  points  between  the  two  data  sets, 
despite  the  relatively  small  number  of  stars  in  common. 

Stars  in  common  are  found  in  four  isolated  regions  that  coin¬ 
cide  with  the  historical  and  well-known  Kapteyn  selected  areas 
113, 92,  95,  and  113.  We  determine  the  zero-point  offsets  between 
the  SDSS  and  Stetson’s  photometry  for  each  region  separately  by 
synthesizing  B  VRI  magnitudes  from  SDSS  gri  photometry  and 
comparing  them  to  Stetson’s  measurements.  The  implied  zero- 
point  errors  (which,  of  course,  can  be  due  to  either  the  SDSS  or  the 
Stetson  data  set,  or  both)  are  listed  in  Table  8.  For  regions  1-3  the 
implied  errors  are  only  a  few  millimagnitudes  (except  for  the  B  —  g 
color  in  region  1 ).  The  discrepancies  are  much  larger  for  the  three 
red  colors  in  region  4.  A  comparison  with  the  results  of  internal 
SDSS  tests  described  in  §§  2.7.1  and  2.7.2  suggests  that  these  dis¬ 
crepancies  are  more  likely  due  to  zero-point  offsets  in  Stetson’s 
photometry  for  this  particular  region  than  to  problems  with  SDSS 
photometry.  We  contacted  P.  Stetson,  who  confirmed  that  his  ob¬ 
serving  logs  were  consistent  with  this  conclusion.  Only  a  small 
fraction  of  stars  from  Stetson’s  list  are  found  in  this  region. 

Given  the  results  presented  in  this  section,  we  conclude44  that 
the  rms  for  the  spatial  variation  of  zero  points  in  the  SDSS  stripe 
82  catalog  is  below  0.01  mag  in  the  gri  bands. 

3.  THE  UTILITY  OF  THE  SDSS  STRIPE  82 
STANDARD  STAR  CATALOG 

As  examples  of  the  use  of  the  standard  star  catalog,  we  discuss 
the  calibration  of  data  obtained  in  nonphotometric  conditions  and 


44  Here  we  assumed  that  it  is  a  priori  unlikely  that  the  SDSS  and  Stetson’s  zero- 
point  errors  are  spatially  correlated. 


a  detailed  and  robust  measurement  of  the  morphology  of  the  stel¬ 
lar  locus  in  color-color  space. 

3.1.  Calibration  of  Nonphotometric  Data 

The  existence  of  a  technique  to  photometrically  calibrate  non¬ 
photometric  data  would  greatly  increase  the  efficiency  of  tele¬ 
scopes.  As  one  particular  example  of  how  our  catalog  can  support 
a  large  project,  consider  the  SDSS-II  supernova  survey  ( Sako  et  al. 
2005 ).  This  survey  aims  to  obtain  repeat  images  of  stripe  82  with 
a  sufficient  cadence  to  enable  the  discovery  of  new  Type  la  super¬ 
novae.  This  requirement  sometimes  results  in  observations  ob¬ 
tained  through  clouds  with  several  magnitudes  of  extinction.  In 
such  highly  nonphotometric  conditions  the  standard  photometric 
calibration  described  in  §  2.2  fails  because  the  fields  with  stan¬ 
dard  stars  are  too  sparsely  distributed  to  be  able  to  resolve  fast 
variations  in  cloud  extinction. 

3.1.1.  A  Method  to  Track  Fast  Cloud  Extinction  Variations 

Due  to  its  high  stellar  density,  the  standard  star  catalog  de¬ 
scribed  in  this  paper  can  be  used  for  calibration  of  data  obtained 
in  grossly  nonphotometric  conditions.  The  typical  number  of  cal¬ 
ibration  stars  in  each  SDSS  field  (9x13  arcmin2 )  at  high  Galactic 
latitudes  is  10-15  in  the  u  band,  40-50  in  the  gri  bands,  and  30- 
40  in  the  z  band.  Based  on  tests  of  several  nonphotometric  SDSS-II 
runs,  it  was  found  that  the  cloud  extinction  variations  can  be  tracked 
with  a  sufficient  temporal  resolution  (~3  s)  to  obtain  photometric 
zero-point  accuracy  comparable  to  that  characteristic  for  photo¬ 
metric  nights  ( 1  %— 2%  in  gri  and  2%— 3%  in  u  andz;  Ivezic  et  al. 
2004a). 


TABLE  8 

Comparison  with  Stetson’s  Standards:  Photometric  Zero-Point  Variations 


Color  ResidualsR1a  crR1b  iVR1c  Residuals^3  crR2b  tVR2°  ResidualsR3a  crR3b  (VR3C  ResidualsR4a  crR4b  NR4C 


B-g .  -29  21  92  6  27  165  8  42  155  -4  27  281 

V-g .  0  17  99  0  15  217  6  25  161  17  19  282 

R-r .  -6  16  58  4  16  135  -8  12  11  39  27  60 

I-i .  -11  16  94  6  18  205  2  16  124  19  15  47 


a  The  median  value  of  residuals  ( in  millimagnitudes)  for  transformations  listed  in  the  first  column,  evaluated  separately  for  regions  1—4,  defined  as  follows:  R 1 ,  R.  A.  ~  325, 
deck  <  0;  R2,  R.A.  ~  15;  R3,  R.A.  ~  55;  R4,  R.A.  ~  325,  deck  >0. 
b  The  rms  scatter  for  the  transformation  residuals  (in  millimagnitudes). 
c  The  number  of  stars  in  each  region  with  good  photometry  in  the  required  bands. 
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The  calibration  is  done  in  two  steps.  First,  the  implied  zero 
points  (whose  variation  is  dominated  by  cloud  extinction),  zp, 
defined  by 


^calibrated  =  -2.5  log  (counts)  +  zp,  (11) 

are  computed  for  each  star  and  median-filtered  in  time  using  a 
window  with  five  stars  in  order  to  avoid  outliers.  Note  that  we 
assume  that  the  clouds  are  gray  and  do  not  allow  for  color  terms, 
an  assumption  which  is  justified  a  posteriori  (see  §3.1 .4).  In  the 
second  step,  zero  points  are  evaluated  for  each  2048  pixel  wide 
(cross-scan  direction)  and  100  pixel  long  (in-scan  direction)  im¬ 
age  segment,  hereafter  called  a  “calibration  patch”  (not  to  be  con¬ 
fused  with  secondary  star  patches,  discussed  in  §  2.2).  That  is,  a 
calibration  patch  is  an  ~9  arcmin2  large  rectangle  with  an  as¬ 
pect  ratio  of  1  : 20,  and  the  zero  points  are  evaluated  every  2.6  s 
( but  note  that  the  variations  are  smoothed  out  by  the  54  s  long 
exposure  time). 

The  patch  is  much  narrower  in  the  in-scan  direction  because 
tests  have  shown  that  zero-point  gradients  across  a  field  are  much 
larger,  by  a  factor  of  10-50,  in  this  direction  (see  Fig.  14).  Con¬ 
sider  three  stars:  star  A;  star  B,  which  is,  say,  25'  (the  column-to- 
column  separation)  away  from  star  A  in  the  scan  direction;  and 
star  C,  which  is  25'  away  from  star  A  in  the  cross-scan  direction. 
Stars  A  and  C  are  observed  at  the  same  time,  and  the  difference  in 
their  implied  zero  points  measures  the  structure  function  (SF)  of 
cloud  opacity  on  a  25'  spatial  scale.  This  is  true  irrespective  of  the 
cloud  motion  relative  to  the  boresight.  Here  the  SF  of  cloud  opac¬ 
ity  is  defined  as  the  rms  width  of  the  distribution  of  zero-point  dif¬ 
ferences  evaluated  for  pairs  of  points  separated  by  some  distance. 

On  the  other  hand,  stars  A  and  B  are  observed  at  times  that  dif¬ 
fer  by  1 .7  minutes.  If  the  component  of  the  cloud  angular  veloc¬ 
ity  on  the  sky  relative  to  the  boresight  and  parallel  to  the  scanning 
direction  is  w°  minute- 1 ,  the  zero-point  difference  for  stars  A  and 
B  samples  the  cloud  structure  on  spatial  scales  of  25 '(uj/ujs), 
where  uis  =  0.25°  minute- 1  is  the  sidereal  scanning  rate  (for  sim¬ 
plicity,  we  assumed  t which  is  supported  by  the  data).  The 
observed  behavior  of  zero  points,  such  as  that  shown  in  Fig¬ 
ure  14,  implies  wind  velocity  in  the  range45  ui  =  3°-15°  minute-1, 
or  uj/ijJs  ~  12-60.  Hence,  drift  scanning  has  the  unfortunate  prop¬ 
erty  that  the  motion  of  an  inhomogeneous  extinction  screen  with  a 
speed  much  larger  than  the  sidereal  scanning  rate  greatly  magni¬ 
fies  the  effective  zero-point  variations  in  the  scan  direction. 

The  zero  points  for  each  calibration  patch  are  computed  by 
taking  all  the  stars  from  the  patch,  or  finding  the  closest  three  stars 
for  sparsely  populated  patches,  and  adopting  the  median  value  of 
their  zero  points.  This  is  certainly  not  the  only,  nor  perhaps  the 
ideal,  approach  to  calibrate  patches,  but  we  found  that  it  works 
well  in  practice.  The  zero-point  error  is  evaluated  from  the  rms 
scatter  of  zp  evaluated  for  each  calibration  star,  divided  by  the 
square  root  of  the  number  of  stars.  We  now  discuss  the  perfor¬ 
mance  of  this  method. 

3.1.2.  Performance  and  Quality  Tests 

The  top  panel  in  Figure  15  summarizes  the  behavior  of  cloud 
extinction  in  the  r  band,  as  measured  by  the  zero  point  zp  dis¬ 
cussed  above,  for  an  SDSS-II  supernova  run  (5646)  obtained  in 
strongly  nonphotometric  conditions.  Although  the  cloud  extinc¬ 
tion  during  the  first  90  minutes  (corresponding  to  150  SDSS  fields) 


45  This  range  is  equivalent  to  angular  speeds  of  up  to  half  of  the  Moon's  diam- 
eter  per  second.  The  plausibility  of  this  wind  velocity  range  was  verified  in  exten¬ 
sive  visual  observations  of  the  full  Moon  during  frequent  grossly  nonphotometric 
nights  in  Seattle. 
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Fig.  14. — Comparison  of  cloud  extinction  gradients  in  the  in- scan  (right  ascen¬ 
sion,  horizontal  axis)  and  cross-scan  (declination,  vertical  axis)  directions  for  SDSS 
run  5759,  on  a  spatial  scale  of  ~0.4°  (each  point  is  derived  using  zero  points  from 
three  calibration  patches).  Note  the  different  axis  scales.  For  this  particular  mn,  the 
distribution  width  is  9.6  times  larger  for  the  in-scan  than  for  the  cross-scan  direc¬ 
tion.  This  is  a  consequence  of  cloud  motion  relative  to  the  boresight  and  the  drift¬ 
scanning  technique.  [See  the  electronic  edition  of  the  Journal  for  a  color  version 
of  this  figure.] 


varies  between  0  and  ~6  mag,  it  is  possible  to  robustly  calibrate 
these  data.  Figure  1 6  zooms  in  on  an  8  minute  stretch  of  the  same 
data  where  the  cloud  extinction  varies  between  0  and  ~3  mag, 
with  changes  as  fast  as  0.05  mag  s-1  (almost  2  mag  per  SDSS 
field).  As  shown  in  the  figure  {third  row,  left  panel),  the  residuals 
have  a  distribution  width  of  only  0.07  mag.  The  right  panels  in 
the  third  rows  ofFigures  15and  16  demonstrate  that  most  of  this 
scatter  is  contributed  by  random  photometric  errors  (i.e.,  errors  in 
extracted  source  counts),  rather  than  by  calibration  errors  ( large 
cloud  extinction  results  in  a  smaller  number  of  calibration  stars, 
as  well  as  in  a  lower  signal-to-noise  ratio  for  those  calibration  stars 
that  are  detected).  Even  with  such  a  large  and  rapidly  varying  cloud 
extinction,  the  zero-point  errors  are  smaller  than  0.05  mag,  with  a 
median  value  of  less  than  0.02  mag.  An  example  of  a  run  with 
somewhat  thinner  and  much  more  stable  cloud  cover  is  shown  in 
Figure  17. 

The  calibration  performance  in  other  bands  is  similar.  For  ex¬ 
ample,  although  the  number  of  calibration  stars  is  smaller  in  the  u 
band  than  in  the  r  band,  the  median  zero-point  error  for  the  same 
stretch  of  data  as  shown  in  Figure  17  is  still  only  0.01  mag,  as  il¬ 
lustrated  in  Figure  18. 

3.1.3.  The  Summary  of  Calibration  Accuracy 

A  summary  of  the  final  zero-point  errors  as  a  function  of 
cloud  extinction  and  band  for  one  of  the  worst  runs  is  shown  in 
the  left  column  of  Figure  19.  As  the  figure  shows,  the  data  can 
be  calibrated  with  small  zero-point  errors  even  for  such  a  bad 
case.  Typically,  the  zero-point  errors  for  the  same  cloud  extinc¬ 
tion  are  about  twice  as  small  as  in  this  run.  A  calibration  sum¬ 
mary  for  a  run  with  optically  thick  but  exceptionally  smooth  clouds 
is  shown  in  the  right  column  of  Figure  19.  Overall,  for  cloud  ex¬ 
tinction  of  X mag,  the  zero-point  uncertainty  is  typically  smaller 
than  (0.02-0. 05)X  for  95%  of the  calibration  patches,  with  a  me¬ 
dian  o/(0.01-0.02)X. 
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Fig.  15. — Top  panel :  Summary  of  the  behavior  of  cloud  extinction  in  the  r  band  over  1 .5  hr  during  SDSS-II  run  5646.  Individual  calibration  stars  are  shown  by  dots, 
and  the  adopted  zero  point  is  shown  by  the  line.  Second  panel.  Calibration  residuals  for  each  star  (dots) .  The  rms  scatter  for  these  residuals  evaluated  for  each  field  is  shown 
by  the  line.  Third  row  of  panels :  Distribution  of  the  residuals  is  shown  in  the  left  panel  (solid  line).  The  median  and  equivalent  Gaussian  a  evaluated  from  the  interquartile 
range  are  also  shown  in  this  panel,  as  well  as  a  Gaussian  corresponding  to  these  parameters  (dashed  line).  The  right  panel  is  analogous,  except  that  the  residuals  are  nor¬ 
malized  by  the  expected  errors.  Bottom  panels :  Distribution  of  implied  cloud  extinction  is  shown  in  the  left  panel,  and  the  distribution  of  standard  errors  for  the  adopted 
photometric  zero  points  (computed  from  the  rms  width  of  the  distribution  of  residuals)  is  shown  in  the  right  panel  (a  few  points,  about  4  out  of  1 800  calibration  patches,  are 
off  scale).  [See  the  electronic  edition  of  the  Journal  for  a  color  version  of  this  figure .] 


3.1.4.  The  Cloud  Color  and  Structure  Function 

We  detect  no  dependence  of  the  calibration  residuals  on  the  stel¬ 
lar  color  or  cloud  thickness  at  the  level  of  a  few  millimagnitudes. 
This  is  consistent  with  the  lack  of  selective  extinction  by  clouds. 
The  lack  of  a  color  correlation  in  the  u  and  z  bands  implies  that  the 
well-known  cloud  grayness  extends  beyond  visual  wavelengths. 
Another  method  to  quantify  selective  extinction  by  clouds  is  to  di¬ 
rectly  compare  zero  points  from  different  bands.  As  shown  in  Fig¬ 
ure  20,  the  cloud  extinction  is  similar  in  all  bands  for  most  fields. 
A  few  cases  where  there  are  deviations  of  a  few  tenths  of  a  mag¬ 
nitude  can  be  easily  understood  as  due  to  temporal  changes  in  the 


cloud  opacity  (recall  that  the  data  from  different  bands  are  ob¬ 
tained  over  ~5  minutes  of  time). 

The  calibration  accuracy  is  determined  by  the  size  of  the  cal¬ 
ibration  patches.  For  example,  a  smaller  patch  would  suffer  less 
from  the  spatial  variation  of  cloud  extinction,  but  it  would  not 
have  enough  stars  to  beat  down  the  noise  of  their  individual  pho¬ 
tometric  measurements  (~0.02  mag  for  sufficiently  bright  stars). 
The  detailed  scaling  of  this  accuracy  with  patch  size  depends  on 
the  cloud  spatial  SF.  The  geometry  of  the  SDSS  camera  allows  us 
to  study  the  cloud  SF  on  scales  exceeding  2°.  Figure  21  compares 
zero  points  in  different  columns  for  two  runs  with  significantly  dif¬ 
ferent  cloud  behavior.  While  zero  points  from  different  columns 
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Fig.  16. — Same  as  Fig.  15,  but  showing  only  ~8  minutes  of  data  with  large  cloud  extinction.  Note  that  the  changes  in  cloud  extinction  are  resolved  down  to  timescales 
well  below  1  minute.  [See  the  electronic  edition  of  the  Journal  for  a  color  version  of  this  figure .] 


generally  track  each  other,  there  can  be  differences  exceeding  a 
magnitude  (they  generally  scale  with  the  cloud  optical  thickness). 
These  differences  increase  with  the  distance  between  the  camera 
columns.  Figure  22  shows  typical  behavior:  for  small  spatial  scales 
(<2°)  the  SF  is  roughly  a  linear  function  of  distance,  and  it  scales 
roughly  linearly  with  the  cloud  extinction.  At  a  1°  scale,  the  SF 
is  typically  of order  2%—10%  of  the  cloud  extinction.  For  exam¬ 
ple,  even  for  clouds  3  mag  thick,  the  SF  at  2'  scales  is  typically 
<0.01  mag. 

3.1.5.  Implications  for  Surveys  Such  as  LSST 

The  LSST  is  a  proposed  imaging  survey  that  will  attempt  to 
maximize  its  observing  time  by  accepting  nonphotometric  condi¬ 


tions.  At  the  same  time,  it  has  adopted  exquisite  requirements  for 
its  photometric  accuracy,  including  1%  accuracy  of  its  internal 
photometric  zero-point  errors  across  the  sky.  Our  analysis  allows 
us  to  answer  the  following  question:  What  is  the  largest  cloud 
extinction  that  still  allows  photometric  calibration  accurate  to 
1%  ? 

A  similar  approach  to  the  calibration  of  LSST  data  as  presented 
here  (assuming  that  a  standard  star  catalog  is  available,  e.g.,  from 
prior  demonstrably  photometric  nights)  would  benefit  from  sev¬ 
eral  effects: 

1 .  The  LSST  will  not  use  a  drift-scanning  technique,  and  thus 
the  calibration  patches  can  be  squares;  for  the  same  area,  this  re¬ 
sults  in  an  ~5  times  smaller  angular  scale  (^3'),  compared  to  the 
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Fig.  1 7. — Same  as  Fig.  15,  but  showing  1 .5  hr  of  data  from  run  5759,  which  had  somewhat  thinner  and  much  more  stable  cloud  cover.  Note  that  the  median  photometric 
zero-point  error  is  below  0.01  mag,  although  the  median  cloud  extinction  is  larger  than  1  mag.  [See  the  electronic  edition  of  the  Journal for  a  color  version  of  this  figure.'] 


1:20  rectangles  we  have  used  to  calibrate  SDSS  drift-scanning 
data.  On  these  angular  scales,  the  cloud  SF  is  roughly  linear,  and 
thus  the  zero-point  error  is  ~5  times  smaller,  or  of  the  order  of 
1%  or  less  through  clouds  as  thick  as  1  mag  (conservatively  as¬ 
suming  that  SDSS  errors  would  be  0.05W;  see  §  3.1.3).  We  note 
that  the  shorter  exposure  time  for  LSST  (30  s,  or  about  a  factor  of 

2  shorter  than  for  SDSS)  is  not  relevant  because  clouds  would 
typically  move  by  more  than  a  degree  during  the  exposure.  This 
is  a  distance  more  than  an  order  of  magnitude  larger  than  the  size 
of  the  calibration  patch,  and  thus  the  SF  analysis  remains  valid. 

2.  LSST  data  will  be  deeper  than  that  of  SDSS  by  about  2— 

3  mag.  With  a  conservative  assumption  that  log  N  oc  0.3  mag  for 
faint  stars  (0.6  for  Euclidean  counts),  the  surface  density  of  the 
calibration  stars  will  be  about  10  times  larger  for  LSST  than  for 


SDSS.  This  larger  density  enables  10  times  smaller  patches,  or 
about  3  times  smaller  angular  scale  for  calibration  (~  1'),  resulting 
in  another  factor  of  3  improvement  of  accuracy. 

3 .  F itting  a  smooth  function  for  cloud  opacity  over  several  cal¬ 
ibration  patches  would  result  in  further  improvements. 

The  first  two  points  predict  thatTXST data  could  be  calibrated 
with  the  required  1%  accuracy  even  through  3  mag  thick  clouds. 
Given  the  various  extrapolations,  we  conservatively  suggest  the 
range  of  1-3  mag  as  the  upper  limit  on  the  acceptable  cloud 
opacity.46 


46  Of  course,  cloud  opacity  decreases  the  imaging  depth,  and  data  with  clouds 
thicker  than  ~  1  mag  may  be  undesirable  for  reasons  other  than  calibration  accuracy. 
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Fig.  1 8 . — Same  as  Fig.  1 7,  but  showing  the  w-band  calibration  summary.  Despite  the  smaller  number  of  calibration  stars  than  in  redder  bands,  and  over  a  magnitude  of  cloud 
extinction,  it  is  still  possible  to  photometrically  calibrate  these  data  with  a  median  error  of  only  ~0.01  mag.  [See  the  electronic  edition  of  the  Journal for  a  color  version  of  this 
figure.) 


3.2.  The  Morphology  of  the  Stellar  Locus 

The  improved  accuracy  of  averaged  photometry  provides 
“crisper”  color-color  diagrams  and  also  reveals  new  morphol¬ 
ogical  features.  An  example  of  such  a  color-color  diagram  is 
shown  in  Figure  23. 

This  is  a  similarplot  to  that  ofFigure  1  of  Smolcic  et  al.  (2004), 
except  that  only  nonvariable  point  sources  are  shown  (note  the 
absence  of  quasars)  and  averaged  photometry  is  used.  The  white 
dwarf/M  dwarf  “bridge”  discussed  by  Smolcic  et  al.  is  clearly 
visible,  as  well  as  the  locus  of  probable  solar  metallicity  giants 


(this  identification  is  based  on  models;  e.g.,  Kurucz  1979),  which 
deviates  from  the  main  locus  at  u  —  g  ~  2.5  and  g  —  r  ~  1 .  Note 
also  the  well-defined  blue  horizontal-branch  locus  ( u  —  g  ~  1.1 
and  g  —  r  from  —0.3  to  0.1)  and  the  white  dwarf  locus  (u  —  g  ~ 
0.35  and  g  —  r  from  —0.3  to  ~0.0).  A  new  locus-like  feature  that 
is  not  visible  in  Figure  1  of  Smolcic  et  al.  is  discernible  at  u  — 
g  ~  0  and  g  —  r  ~  —0.2.  The  great  value  of  the  accurate  M-band 
photometry  is  clearly  evident  at,  e.g.,  g  —  r  =  —0.2;  the  u  —  g 
color  distribution  is  trimodal!  The  bluest  branch  is  consistent 
with  He  white  dwarfs  and  the  middle  branch  with  hydrogen  white 
dwarfs,  as  supported  by  Bergeron  et  al.  (1995)  white  dwarf  models 
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Fig.  19. — Summary  of  calibration  accuracy  as  a  function  of  cloud  extinction  and  band  ( urz ,  as  marked  in  the  panels).  Each  small  symbol  represents  one  calibration 
patch  (an  ~9  arcmin2  large  rectangle  with  a  1:20  aspect  ratio).  The  zero-point  error  is  determined  from  the  rms  of  the  photometric  residuals.  The  large  symbols  show  the 
median  zero-point  error  in  1  mag  wide  bins  of  cloud  extinction.  The  left  column  shows  data  for  one  of  the  photometrically  worst  SDSS-II  supernova  runs  (5646),  and  the 
right  column  shows  data  for  a  run  with  optically  thick  but  exceptionally  smooth  clouds  (5759).  Note  that  the  data  can  be  calibrated  with  zero-point  errors  typically  smaller 
than  a  few  percent  even  through  clouds  several  magnitudes  thick.  [See  the  electronic  edition  of  the  Journal  for  a  color  version  of  this  figure .] 


and  detailed  analysis  of  SDSS  spectra  ( Eisenstein  et  al.  2006). 
The  reddest  branch  is  made  of  blue  horizontal-branch  stars  (see 
Sirko  et  al.  2004  and  references  therein). 

The  exciting  fact  that  one  can  distinguish  He  and  H  white 
dwarfs  using  photometry  alone  is  a  consequence  of  the  improved 
photometric  accuracy  due  to  averaging  many  epochs.  Figure  24 
reiterates  that  point.  Note  the  striking  difference  between  the  two 
bottom  panels:  while  one  could  be  convinced  that  the  He  white 
dwarf  sequence  is  a  real  feature  in  the  bottom  right  panel,  its  ex¬ 
istence  is  clearly  evident  when  using  the  improved  photometry, 
as  shown  in  the  bottom  left  panel.  In  summary,  the  multiepoch 
observations  provide  both  the  identification  of  variable  sources 
and  much  more  accurate  colors  for  nonvariable  sources.  This 
bodes  well  for  science  deliverables  from  upcoming  large-scale 
imaging  surveys.  For  example,  LSST  will  obtain  over  its  10  yr 
mission  repeat  imaging  similar  to  that  discussed  here,  but  about 
2.5  mag  deeper,  with  about  100  or  more  observations  per  band 


and  object  and  over  an  area  about  2  orders  of  magnitude  larger. 
Although  these  new  surveys  will  not  have  a  spectroscopic  com¬ 
ponent  like  SDSS  did,  the  multiepoch  nature  of  their  imaging 
will  provide  alternative  information-rich  data  sets. 

4.  DISCUSSION  AND  CONCLUSIONS 

Using  repeated  SDSS  measurements,  we  have  constructed  a 
catalog  of  over  1  million  candidate  standard  stars.  The  catalog  is 
publicly  available  from  the  SDSS  Web  site.47  Several  indepen¬ 
dent  tests  suggest  that  both  internal  zero-point  errors  and  random 
photometric  errors  for  stars  brighter  than  19.5,  20.5,  20.5,  20, 
and  18.5  in  ugriz,  respectively,  are  at  or  below  0.01  mag  (about 
2-3  times  better  than  for  an  individual  SDSS  scan).  This  is  by  far 
the  largest  existing  catalog  with  multiband  optical  photometry 


47  See  http://www.sdss.org/dr5/products/value_added/index.html. 
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Fig.  20. — Color  of  cloud  extinction  in  SDSS  bands.  Each  symbol  represents  one  field  and  shows  the  difference  in  cloud  extinction  between  the  two  bands  as  a  function 
of  the  r-band  extinction.  The  measurements  in  different  bands  are  obtained  over  ~5  minutes  of  time,  and  thus  even  gray  clouds  with  spatially  varying  extinction  could  pro¬ 
duce  the  observed  nongray  (nonzero)  values.  The  dashed  lines  in  the  bottom  two  panels  indicate  the  expected  correlation  if  the  color  variations  are  due  to  temporal  changes 
in  the  gray  cloud  thickness  (rather  than  due  to  intrinsic  color  changes).  [See  the  electronic  edition  of  the  Journal  for  a  color  version  of  this  figure .] 
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Fig.  21. — Comparison  of  cloud  extinction  independently  measured  for  six 
camera  columns.  Top  two  panels :  30  minutes  of  measurements  in  the  g  and  r  bands 
for  the  same  run  (5759)  shown  in  Fig.  17.  Individual  camera  columns  are  color- 
coded  according  to  the  legend  shown  in  Fig.  6  {top).  Third  panel :  Difference 
between  the  r-band  zero  points  measured  in  one  of  the  edge  columns  (6 )  and  zero 
points  from  the  other  five  columns.  Bottom  two  panels :  Same  as  second  and  third 
panels,  but  with  data  from  a  run  with  exceptionally  patchy  clouds  (5646;  the  first 
30  minutes  of  the  data  from  Fig.  1 5  are  shown).  Note  the  varying  scale  for  they- axis. 

accurate  to  ~  1  %  and  breaks  the  accuracy  barrier  discussed  by, 
e.g.,  Stubbs  &  Tonry  (2006  and  references  therein;  see  also  Maiz 
Apellaniz  [2006]  for  a  discussion  of  1  %  photometry  in  other  pho¬ 
tometric  systems).  These  observations  were  not  obtained  for  the 
specific  purpose  of  calibration  but  were  part  of  the  regular  SDSS 
observational  program.  When  compared  to,  for  example,  the  he¬ 
roic  calibration  efforts  by  A.  U.  Landolt,  P.  Stetson,  and  others,  it 
seems  justified  to  call  the  method  presented  here  “industrial” 
photometry.  However,  the  catalog  presented  here  is  not  without 
its  problems. 

The  selection  of  candidate  stars  was  simply  based  on  the  ab¬ 
sence  of  variability.  It  is  fairly  certain  that  not  all  variable  sources 
are  recognized  because  of  the  limited  number  of  repeated  ob¬ 
servations  (~10).  For  example,  an  eclipsing  binary  with  a  much 
shorter  eclipse  duration  than  the  orbital  period  could  easily  es¬ 
cape  detection.  Furthermore,  some  of  these  sources  may  not  even 
be  stars.  A  cross-correlation  with  the  SDSS  spectroscopic  data¬ 
base  yields  99,000  matches  in  the  overlapping  region,  including 
44,000  spectra  classified  as  stars.  About  70  candidate  standard 
stars  are  actually  spectroscopically  confirmed  quasars!  Appar¬ 
ently,  a  small  fraction  of  quasars  (a  few  percent;  for  a  detailed 
analysis,  see  B.  Sesar  et  al.  2007,  in  preparation)  cannot  be  de¬ 
tected  by  variability  (at  least  not  using  the  number  of  epochs,  their 
time  distribution,  and  the  photometric  accuracy  employed  in  this 
work).  Indeed,  we  have  also  found  three  spectroscopically  con¬ 
firmed  SDSS  quasars  among  Stetson’s  standards  which  were  ob¬ 
served  20-30  times  and  which  showed  no  variation.  Similarly, 
about  300  candidate  standard  stars  have  SDSS  spectra  classified 
as  galaxies.  Nevertheless,  the  inspection  of  color-color  diagrams 


distance  (deg) 

Fig.  22. — Cloud  SF  in  the  r  band  for  run  5759  and  for  the  same  stretch  of  data 
as  shown  in  Fig.  17  (the  median  cloud  extinction  is  1 .3  mag).  The  circles  show  the 
rms  width  of  the  distribution  of  zero-point  differences  between  camera  column  1 
and  the  five  other  columns.  This  width  is  corrected  for  a  0.015  mag  contribution 
from  the  measurement  errors  and  shown  by  squares.  The  triangles  show  the  width 
of  the  distribution  of  zero-point  differences  in  the  in-scan  direction,  with  the  dis¬ 
tance  scale  multiplied  by  30.  This  multiplication  factor  measures  the  cloud  speed 
relative  to  the  boresight  in  the  in-scan  direction  (see  §  3.1.1).  [See  the  electronic 
edition  of  the  Journal  for  a  color  version  of  this  figure .] 

strongly  suggests  that  the  overwhelming  majority  of  the  standard 
stars  are  found  on  the  stellar  locus. 

Remaining  systematic  errors  are  another  important  concern. 
Effectively,  we  have  assumed  that  PT  problems  average  out  in 
many  patches  when  deriving  flat-field  corrections  using  stellar 
colors.  This  may  not  be  true  at  a  level  not  much  smaller  than  1%, 
and  thus  the  remaining  gray  problems  at  such  a  level  may  be  pres¬ 
ent  in  the  catalog.  Despite  these  residual  problems,  we  believe  that 


u-g 


Fig.  23. — The  g  —  r  vs.  u  —  g  color-color  diagrams  for  all  nonvariable  point 
sources  constructed  with  the  improved  averaged  photometry  {dots).  Various  stellar 
models  (Kurucz  1979;  Bergeron  et  al.  1995;  Smolcic  et  al.  2004)  are  shown  by 
lines,  as  indicated  in  the  figure. 
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Nobs>3,  non-variable  Single-epoch  observations 


Fig.  24. — Advantages  of  repeated  photometric  measurements.  Top  left  panel:  Blue  comer  of  the  g  —  r  vs.  w  —  g  diagram  from  Fig.  23  for  all  point  sources  with  the  aver¬ 
aged  photometry.  Top  right  panel :  Only  the  variable  sources  (dominated  by  low-redshift  quasars).  Bottom  left  panel:  Nonvariable  sources  (dominated  by  white  dwarfs) 
classified  using  low-order  light-curve  moments.  Bottom  right  panel:  Same  nonvariable  sources,  but  using  their  DR5  single-epoch  photometry.  A  comparison  of  the  bottom 
two  panels  shows  the  striking  improvement  made  possible  by  the  use  of  multiple  observations  of  the  same  field. 


internal  consistency  of  the  catalog  (i.e.,  when  ignoring  A„,  from 
eq.  [4])  is  such  that  the  rms  width  for  the  function  6„,(R.  A. ,  decl.) 
from  equation  (4)  evaluated  for  all  stars  in  the  catalog  is  at  most 
0.01  mag  in  the  gri:  bands  and  perhaps  just  slightly  larger  in  the  u 
band  (very  unlikely  exceeding  0.02  mag).  In  addition  to  gray  prob¬ 
lems  and  overall  flat-field  errors,  the  dependence  of  flat  fields  on 
source  color  is  probably  the  largest  remaining  systematic  error. 

We  illustrate  several  uses  for  this  catalog,  including  the  calibra¬ 
tion  of  highly  nonphotometric  data  and  robust  selection  of  stars 
with  peculiar  colors.  We  find  that  LSST  and  similar  surveys  will 
be  able  to  observe  in  partially  cloudy  (nonphotometric)  nights  be¬ 
cause  even  cloudy  data  can  be  accurately  calibrated  with  a  suffi¬ 
ciently  dense  network  of  calibration  stars.  Such  a  dense  network 
will  be  self-calibrated  by  LSST  very  soon  after  first  light,  using  an 
approach  developed  for  SDSS  data  by  Padmanabhan  et  al.  (2007). 
Given  such  a  network,  SDSS  experience  suggests  that  LSST  can 
maintain  its  required  photometric  calibration  accuracy  of  1  %  even 
when  observing  through  1-3  mag  thick  clouds. 

Perhaps  the  most  exciting  conclusion  of  this  work  is  that  it  may 
become  obsolete  in  only  a  few  years  due  to  the  advent  of  next- 
generation  surveys  such  as  Pan-STARRS  and  LSST. 
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